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INTRODUCTION 


This  Annual  Progress  Report  contains  investigations  on  detonation 
driven  magnet ogasdynamic  and  solid  state  power  generators  and  related 
research  problems.  The  work  was  conducted  under  ONR  Contract  N00014- 
81-WR- 10107  in  the  period  from  February  1  to  September  30,  1981. 

The  electric  currents  and  voltages  produced  by  electromagnetic 
induction  in  explosion  driven  magnetogasdynamic  generators  with  homo¬ 
geneous  external  magnetic  fields  Bq  are  calculated  for  two  shock  flow 
models,  (1)  the  (implosion  produced)  jet  flow  v(t)  *  voH(t)  with  con¬ 
stant  shock  speed  vq,  and  (2)  the  plane  detonation  flow  with  shock 
-1/3 

speed  v(t)  «  t  (similarity  solution).  External  load  circuits  with 


resistance  R  and  inductance  L  are  considered  which  are  connected 
o  o 

i)  to  the  upstream  and  ii)  to  the  downstream  ends  of  the  generator 
electrodes  so  that  the  magnetic  self-fields  B  of  the  electric  currents 


are  i)  antiparallel  and  ii)  parallel  to  the  external  tic  field 

It  is  shown  that  optimum  power  is  generated  for  generators  with 


jet  shock  flow  and  load  circuits  of  type  (ii),  i.e.,  positive  super¬ 


position  of  induced  and  external  magnetic  fields.  In  this  case,  elec-  X* 

5  /4%  *  Vn 

trie  energies  AE  *  10  Joule  can  be  generated  with  external  fields  X  X  X  > 


Bq  »  1  Tesla  per  detonation.  The  importance  of  the  results  in  the  ^ 


COP* 


XX 

/  '*e 


\V\ 

v  . 


Chapters  II  and  III  for  the  design  of  detonation  driven  generators  is 
obvious . 

One  of  the  main  technical  problems  of  explosion  driven  magneto- 
g  a  adynamic  generators  is  the  requirement  of  large  external  magnetic 
fields  Bq  >  1  Tesla,  which  are  difficult  to  procure  in  connection  with 
military  applications.  For  this  reason,  a  solid  state  power  generator 
is  analyzed  in  which  the  driving  electromotoric  force  is  due  to  electric 
polarization  in  the  stress  field  of  an  explosion  produced  shock  wave# 

The  solid  is  assumed  to  be  a  one-dimensional  slab  between  plane  elecr- 
trodes,  which  are  connected  through  an  external  load  circuit •  The  tem¬ 
poral  current  build-up  in  the  external  circuit,  due  to  electric  polari¬ 
zation  and  electric  conduction  behind  the  shock  front  of  the  stress 
wave,  space  charge  and  electric  polarization  relaxations,  is  determined 
by  an  inhomogeneous  Volterra  integral  equation.  For  shocked  solids  with 
cross  section  of  the  order  of  one  square  meter,  energies  of  the  order 
AE  -  10^  Joule  per  pulse  can  be  achieved.  This  scheme  for  an  electric 
generator,  which  does  not  require  external  magnetic  fields,  appears  to 
be  promising  for  experimental  Investigation  (Chapter  IV) • 

An  other  method  of  avoiding  large  external  magnetic  fields  Bq  > 

1  Tesla  is  to  build  an  explosion  driven  MOD  generator,  in  which  a  weak 

-2 

external  magnetic  field  Bq  -  10  Tesla  is  compressed  between  the  plasma 

shock  front  (high  magnetic  Reynolds  numbers  R  «  y  ova)  and  a  copper 

m  o 

liner.  The  exact  relativistic  theory  of  the  flux  compression  MGD  gen¬ 
erator  leads  to  two  coupled  integrodlfferentlal  equations  which  appear 


2 


to  be  unsolvable  by  analytical  methods.  For  this  reason,  the  electromagne¬ 
tic  diffusion  and  induction  in  accelerated  conductors  in  external  ^  fields 

o 

was  analyzed  in  the  parabolic  diffusion  approximation,  which  is  valid  for 

—18  7  “1 

small  relaxation  times  x  *  eQ/a  (t  -  10  sec  for  copper,  a  -  6x10  0  /m). 

The  conventional  boundary  conditions  (which  violate  the  continuity  of  the 
tangential  electric  field  )  were  replaced  by  new  boundary  conditions,  which 
consider  the  electromagnetic  transients  outside  of  the  conductor.  Thus,  a 
realistic  evaluation  of  magnetic  flux  compression  and  magnetic  flux  losses 
through  the  plasma  piston  and  the  copper  liner  is  accomplished  (Chapters  V 
and  VI). 

In  Chapter  VII,  a  quantum-kinetic  theory  of  the  anomalous  plasma 
heating  by  high-frequency  electromagnetic  fields  is  presented  based  on  the 
Maxwell-Vlasov  and  Schroedinger  equations.  This  investigation  has  been 
carried  through  in  collaboration  with  visiting  Prof.  S.H.  Kim. 

The  work  on  flux  compression  electric  power  generators  will  be  com¬ 
municated  in  the  1982  Annual  Report.  There,  also  new  types  of  explosion 
driven  generators  will  be  discussed,  which  do  not  require  an  external 
magnetic  field. 
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DETONATION  DRIVEN  INDUCTION  GENERATOR 


WITH  ANTIPARALLEL  EXTERNAL  AND  INDUCED  MAGNETIC  FIELDS 


H.E.  Wilhelm 

Mlchelson  Laboratory,  Naval  Weapons  Center, 
China  Lake,  CA  93555 


Abstract 


Based  on  Maxwell's  equations,  the  electric  current  and  voltage 

pulses  induced  in  explosion  driven  induction  generators  (R  m  u  avL  >>  1) 

m  o 

with  plane  electrodes  and  homogeneous  external  magnetic  field  "Sq  are  cal¬ 
culated  analytically  for  two  plasma  shock  flow  models,  (!)  the  jet  flow 

with  shock  speed  v(t)  -  vQH(t)  and  (ii)  the  plane  detonation  flow  with 

1/3  -1/3 

shock  speed  v(t)  *  (2/3)(Eq/po)  t  .  The  external  load  circuit  with 

resistance  R  and  inductance  L  is  connected  to  the  flow  entrance  ends  of 
o  o 

the  electrodes  so  that  the  magnetic  self-field  B(t)  of  the  generated  cur¬ 
rent  I(t)  is  antiparallel  to  the  transverse  external  field  It  is  shown 

that  the  jet  flow  produces  large  current  and  voltage  pulses  with  extended 
plateaus,  whereas  the  detonation  flow  generates  considerably  smaller  and 
shorter  current  and  voltage  pulses*  In  both  cases,  the  magnetic  self-field 
B(t)  of  the  generated  current  I(t)  is  of  the  order-of-magnltude  of  and 
its  negative  superposition  on  the  external  magnetic  field  is  not  advanta¬ 
geous  for  power  generation. 
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INTRODUCTION 


In  contrast  to  the  conductive  magnetohydrodynamic  generator  with  steady-state 

magnetic  and  flow  fields  ,  induction  generators  are  based  on  the  electromagnetic 

interaction  of  a  conducting  flow  field  with  a  magnetic  field,  where  at  least  one  of 

2) 

the  fields  is  transient  —  so  that  voltages  are  induced  due  to  temporal  flux  changes 
d$/dt  0.  In  some  induction  generators,  a  detonation  plasma  with  a  decreasing  shock 
velocity  v(t)  0  for  t  -►  ®  is  directly  injected  into  the  inter-electrode  space  — \ 
whereas  in  more  recent  experiments^2^ discontinuous  plasma  jets  are  shot  across  an  ex¬ 
ternal,  constant  magnetic  field.  The  plasma  jets  v(t)  -  VQH(t)  of  nearly  constant 
shock  speed  vq  are  produced,  e.g.,  by  implosion  of  a  dense  noble  gas  between  flyer 
plates  and  subsequent  expansion  through  a  diaphragm  into  the  generator  channel  — ^ ,  or 
by  means  of  a  tube  like  arrangement  of  the  solid  explosive  resulting  in  a  cylindrical 
implosion  with  large  axial  shock  velocity  of  the  filling  gas  — \ 

4  5 

Detonation  driven  Induction  generators  have  shock  velocities  v  *  10  -  10  m/sec, 

4  -1 

external  magnetic  fields  B  ■  1  -  10  Tesla,  conductivities  a  >  10  ft  /m,  and  flow 

1  2  10  12  ”3  4  2_  —3 

energies  pv  "  10  -  10  Joule  m  at  pressures  p  >  10  oar  (p  >  10  Kg  flT  ) .  The 

corresponding  magnetic  Reynolds  number  Rffl  is  large,  whereas  the  magnetic  interaction 

number  R^  is  small. 


R_*  U  ovL  »  1 

0  O 


R  -  B2/u  pv2  «  1 
1  o  o 


where  L  -  |V*1/B|.  For  Rn  >>  1»  the  electric  current  Is  restricted  to  a  thin  layer 

6  ~  (t/u  a)*^2  behind  the  shock  front  (t  ■  transit  time  of  flow  in  B  ).  In  view  of 
o  o 

1  2  2  &  f 

the  large  dynamic  pressure  -j  pv  compared  to  the  magnetic  pressure  Bq/2uq  *  10  -  10 

-3 

Joule  a  ,  the  efficiency  of  energy  transformation  is  small,  while  the  performance  is 

A 

significant  since  currents  I  -  (B  /u  )b  >  10  Amps  (b  -  electrode  width)  and  energies 

— *  4  45) 

Ult  >  10  Joule  have  been  produced  experimentally  per  detonation  — ' . 
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Herein,  the  initial -value  problem  for  a  detonation  driven  generator 

is  solved  in  closed  form  by  means  of  Maxwell’s  equations,  in  which  the  flow 

entrance  ends  of  the  coplanar  electrodes  are  connected  by  an  external  L  -R 

o  o 

load  circuit  (Fig.  1)  so  that  the  magnetic  field  B(t)  of  the  generated  cur¬ 
rents  I(t)  is  antiparallel  to  the  external  magnetic  field  Bq.  Two  types  of 
plasma  shock  flows  are  considered  corresponding  to  i)  the  idealized  jet 
flow  with  constant  shock  speed  vq  and  ii)  the  plane  detonation  flow  due  to 
an  explosive  energy  release  Eq  per  unit  area  in  a  gas  of  initial  density  pq: 


v(t)  -  v  H(t)  ,  x(t)  »  v  tH (t )  , 

o  o 

0  <  t  <  t  ,  t  -  x  /v  ; 

00  Q 

and 

v(t)  -  (2/3)(Eo/po)1/3t'1/3  ,  x(t)  -  (Eo/po)1/3t2/3 
0<t<t  ,  t  -  (Eo/po)'1/2x3/2 

where  x(t)  is  the  shock  front  position,  and  t  is  the  transit  time  of  the 
shock  front  x(t)  in  the  duct  of  effective  length  (Fig.  1).  The  self¬ 
similar  solution  (3)  is  asymptotically  correct  for  larger  times  t  ,  but 
can  be  used  for  t  >  0  since  the  singularity  of  v(t)  at  t  *  0  is  without 
physical  consequences  for  generator  applications. 

The  external  circuit  is  assumed  to  be  located  at  x  ■  -x  (Fig.  1), 

o 

where  x  is  an  effective  distance  which  can  be  calculated  once  the  geometries 
o 

of  Che  Inductance  L  and  resistance  R  are  specified. 


(2) 


(3) 


ield  Bq,  and  Rq  -  L  load  circuit  AD.  Shock  front 


INITIAL-VALUE  PROBLEM 


Fig.  1  exhibits  a  one-dimensional  model  for  a  detonation  driven  generator,  in 
which  the  shock  wave  flow  v(t)  *  d5fc(t)/dt  is  initiated  at  time  t  =  0  in  the  plane 
x  ■  0  by  explosive  charges  in  the  space  -xq  <  x  <  0  so  that  the  shock  front  is  at 
x(t)  >  0  for  t  >  0.  The  Cu-electrodes  in  the  planes  z  *  ±c/2  have  the  extensions 
-x  <  x  <  x  <  «  and  -b/2  <  y  <  +b/2,  and  are  connected  by  an  external  circuit  AD 

O  co 

with  load  R  and  inductance  L  .  As  a  result  of  the  v*S  interaction  of  the  conducting 
O  0  o 

shock  front  (a  >  10^0  Vm)  with  the  external  magnetic  field  ^  *  {0fBo,0},  a  current 

I (t)  is  induced  behind  the  shock  front  BC  in  a  cross  section  AxAy  *  b 6  where 
6  -  (t/yQa)  is  the  (small)  width  of  the  current  layer.  The  current  I(t)  flows 
in  the  closed  loop  A-*-B-*-C-»-D-*-A,  where  the  leg  BC  at  5t(t)  moves  with  the  shock 
speed  v(t)  ■  dx(t)/dt. 

Application  of  the  surface  curl  equation  nx[S]  =  y  "t(t)/b  to  the  interfaces  AB, 
BC,  CD,  DA  demonstrates  -  in  the  one-dimensional  approximation  neglecting  end  effects 
at  y  *  ±b/2  and  x  *  -x-  -  that  the  current  I(t)  produces  a  magnetic  self-field  in  the 
extending  loop  volume, 

f(t)  -  {0,-poI(t)/b,0}  ,  -xQ  <  x  <  x(t)  .  |y|  <  j,  M  <  f  t  W 

which  is  antiparallel  to  the  external  magnetic  field  1SQ.  The  total  magnetic  field 
within  the  loop  volume  is 

l(t)  -  {0,Bo-yoI(t)/b,0}  ,  -xQ  <  x  <  x(t)  ,  | y |  <  j  ,  |z|  <  j  ,  (5) 

whereas  $(t)  >1  outside  of  this  region.  Integration  of  the  Maxwell  equation  7x1  - 
o 

-3$/3t  over  the  cross  secti'  ABCD  of  ie  extending  current  loop  yields  by  Stokes' 
theorem 

dr  -  -  +  <f>  <$*%)•& 
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(6) 


where 


E  «  E  +  v*B  ,  v  <<  1/v  e 

o  o 


and  v(t)  is  the  velocity  field  of  the  sections  of  the  closed  loop,  at  the  point 
corresponding  to  the  path  element  dr[v(t)  *  e^dx(t)/dt  for  BC»  but  v(t)  ■  t)  for  CD, 
DA,  AB] .  The  surface  element  dS  *  -e^cdx  forms  with  dr  [pointing  everywhere  in  the 
direction  l(t) ]  a  right  hand  system.  By  Eqs.  (6)  and  (7) 

<»„  +  *,>««  +  f °  Tp*  +  ?'*  1*0  + 

-x  -c/2 

o 

since  /  E *dz  »  IRp  (Rp  "  plasma  resistance),  /  Edz  *  RQI  +  LQdI/dt,  and  /Edz  *  0 

®  D  8  —l 

for  the  sections  AB  and  CD  of  the  highly  conducting  electrodes  (o  ~  10  SI  /m). 


Accordingly, 


where 


(*„  +  VI(t>  +  it { tL0  +  i<c))I(t))  ■ 


L(t)  *  u  clx  +  x(t)]/b  >0  ,  0  <  x(t)  <  x 

o  o 


is  the  time-dependent  inductance  of  the  electrode  sections  AB  &  DC  of  increasing 
length  Ax  -  xo  +  x(t).  Multiplication  of  Eq.  (9)  by  I(t)  leads  to  the  conservation 


<Ro  +  Vl2  +  dT{2  [Lo  +  L(t)11 


2}  +  V  .  t  && 2  b  c 

1  2b  1  dt  dt  oc 


according  to  which  the  ohmic  power  dissipation,  the  temporal  change  of  the  magnetic 

A 

energy,  and  the  work  on  the  current  sheath  per  unit  time  equal  the  power  input  UI. 
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The  plasma  resistance  *  Rp(t)  a  s^ow^y  varying  function  of  t  and  can  be 
calculated  only  via  a  detailed  physical  model  of  the  current  sheath  6  at  x  %  x(t) . 
However,  since  Rq  is  sufficiently  large  in  actual  experiments,  Rq  +  Rp  can  be  treated 
as  a  quasi- con st ant , 

R0  +  Rp(t)  *  constant,  Rq >  Rp  ;  RQ  +  Rp(t)  =  Rq  ,  Rq  »  Rp  .  .  (12) 

For  a  generalized  formulation,  dimensionless  dependent  and  Independent  vari¬ 
ables  are  Introduced  by 


J(t)  *  I(t)/Io  ,  t  =  t/tQ  ,  £(t)  “  x(t)/c 


lo  ’  <Bo/“o)b  •  eo  ■  <“oc2/b)/«o  +  V  ■  ‘o’  (VC>  +  V<V2/b>'<14> 

By  Eqs.  (9)  -  (10)  and  (13)  -  (14),  the  electric  current  I(t)  in  the  explosion  driven 
generator  is  determined  by  the  inhomogeneous  initial-value  problem  with  variable  co¬ 


efficients: 


~  {[Ao  +  C(T)]J(T)}  +  J(x)  - 


J  (t  =  0)  =*  0 


,  0  <  T  <  T 


Eqs.  (15)  -  (16)  are  solved  for  the  jet  and  detonation  flows  in  Eqs.  (2)  and 

(3),  respectively,  which  define  the  dimensionless  coordinate  £(t)  of  the  shock  front. 

12  2 

In  view  of  the  large  flow  energy,  y  pv  >>  Bo/2pq,  it  is  not  necessary  to  correct 
the  above  generator  model  for  velocity  decreases  due  to  Lorentz  and  viscous  forces. 
Furthermore,  viscous  boundary  layers  are  disregarded  since  the  plasma  moves  in  form 
of  a  rectangular  velocity  slug  through  the  duct  in  large  Reynolds  number  flow, 

R  -  pvL/y  >>  1. 
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JET-FLOW  SOLUTION 


For  the  jet  flow,  the  dimensionless  shock  front  coordinate  £(x),  and  the 
transit  time  x  are  by  Eqs.  (2)  and  (13) 

2 

v  y  c  /b  x  R  -HL, 

C(t)  *  CqT  ,  *  c  (R  +r  )  »  T  2,  ^  '  ^ 

O  r  O  UQC  /b 

Since  C(t)  >  0  for  x  >  G,  the  initial-value  problem  (15)  -  (16)  is  without  singu¬ 
larities,  and  has  the  closed-form  solution: 


E0  -(l+£  )/£ 

J(t)  ■  H  -  a  +  x2  T>  ]  ’  o<t<t 

so  o 


The  transit  time  in  Eq.  (17)  satisfies  the  inequality. 


VV* "  T/[~r  +  ~ tt1  >s>  1  •  for  x-  >>  x0  •  Lo 5  (xo/c)lV2/b  ,  (i9) 

V  /b 

which  indicates  that  the  independent  variable  (5o/Aq)x  assume  large  values. 

Eq.  (18)  has  the  limits: 


J(t)  -  (£  /A  )t  »  (£ft/Aft)t  <<  1 
o  o  o  o 


,  (20) 


J(T)  -  £  /( 1  +  £  )  ,  (£  7A  )t  »  1 

o  o  o  o 


(21) 


Accordingly,  the  maximum  current  is  reached  only  asymptotically  for  large  times, 
since  J(x)  is  a  monotouically  increasing  function  of  0  <  x  <  x, 


J  -  Ud  +  £J  ,  I  “  r2  be  /(I  +  £  ),  (£  /A  )t  »  1 

max  o  o  max  o  o  o  o 


ii 


The  generated  voltage,  U  *  RqI  +  LQdI/dt,  is  by  Eq.  (18)  in  dimensionless  represents- 


£  £  -(1+CJ/C.  L  /R  £  £n  -(l+2£  )/£ 

-itt11-  d+rT)  ]  +T-^r(1  +  rT) 


where  U  *  I  R  .  Eq.  (23)  has  the  asymptotic  values 
o  o  o 


L  /R  1  +  2£  £  L  /R  E  £ 

U(x)  *  [1 - 2 — o  ( — - o}]_o  x  +  _2_o  _o  ,  _2.  T  «  i 

o  o  o  o  o  o 


,  (24) 


U(T)  "rrr  *  rT>>1 

o  o 


(25) 


For  1  +  (£  /A  )t  -  (L  /R  t  ) (1  +  2£  )/A  ,  U(r)  assumes  a  minimum, 
o  o  o  o  o  o  o 


E  L/R  1  +  2£  -(l+£  )/£ 

O  r  O  O  /  On  i  O  O 


min  1  +  £ 


-a —  ;  i _ r-u  ( - ] 

+  £  1 X  1  +  2£  1  t  '  A  ' J 


Accordingly,  the  maximum  voltage  value  is  reached  only  asymptotically  for  large  times, 
«o/Ao)T  >>  1’ 


U  -  £  /(I  +  £  )  ,  U  *IR  £  / (1  +  £  ) 

max  o  o  max  o  o  o  o 


Example.  For  b  *  10  3m,  c  ■  1  m,  x  ■  10  ^m,  x^  ■  1  m,  o  ■  1030  ^/m.  Ip  ■ 

o-1c/b6  "  10~^0,  R  +  R_  ■  2xl0_1n,  L  ■  4*«10~®Henry ,  v  *  10 ^m/ sec,  B  *  1  Tesla, 
’op  o  o  o 

-5  -  -4 

the  characteristics  of  the  generator  are:  tQ  *  2irxio  sec,  t  *  10  sec,  £q  ”  2ir/10, 

Aq  -  11/10,  I  -  106/4tt  Amps,  I  •  3xl04  Amps,  1  6xl03Volt,  (01)^  -  2x10®  Watt, 

_ *  L 

AE  ■  Ult  -  10  Joule. 
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Fig,  2  shows  the  current  J(t)  versus  (£  /A  )x  for  £  *  2tt/10  and  £  *  9/10 

o  o  o  o 

based  on  Eq.  (19).  The  distribution  J(x)  rises  steeply  within  0  <  (£o/Aq)x  <  1  and 

then  flattens  out  into  a  shightly  increasing  "plateau"  for  larger  times  x.  In  ex- 
4  5) 

periments  — J ,  the  "plateau"  is  more  constant  or  even  decreases  somewhat  with  in¬ 
creasing  time  t,  which  appears  to  be  due  to  i)  shock  speeds  v(t)  decreasing  slightly 
with  time  t,  and  ii)  external  circuit  leads  to  the  center  of  the  electrodes. 


DETONATION-FLOW  SOLUTION 


The  dimensionless  shock  front  coordinate  £(t),  and  the  transit  time  t  for 
the  detonation  flow  follow  from  Eqs.  (3)  and  (13)  as 


%  .  2/3  .  ,Eo,l/3  -1,  V  /b.2/3  -  „-3/2  3/2 

e(r)  -  c0t  ,  -  (— )  c  (-TTTJ  ’  T  "  5o  (T> 


Since  £(x)  >  0  for  t  >  0,  the  solution  of  the  initial-value  problem  (15)  -  (16)  is 
representable  in  terms  of  an  indefinite  integral  with  regular  integrand: 


J(x)  «  1  - 


.  x  3C*1  [  T1^3-a~^arctan(o()T^^3) ] 

1  *  Ap  i0e _ * 

..  .  2  2/3. 

U  +  aQx  )e 


,  0  <  t  <  t  ,  (29) 


where 


«„  *  W 


Eq.  (29)  is  an  integral  functional  of  aQt  ,  where 


v1/3  -  nT*  — lr1)1/2  ”  1  ■  for  *-  >>  *0  ■  Lo  s  (lto/c)“oc2/b- 

p  c  /b 

°  (31) 

For  small  and  large  times,  Eq.  (29)  gives  the  limits 


t,  ^  2  2/3  1/3  _  . 

J(x)  “  o  r  ,  a  x  <<  1 


,  (32) 


J(t)  -  0 


,  a  x1'3  »  1 
9  o 


J(t)  has  a  true  maximum  J  at  time  t  ■  t  where  dJ(x)/dT  *  0  for  x  *  f .  By  Eqs.  (15), 

max  ^  7 

(28),  and  (29),  J  and  f  are  given  by 


■  "  ~  V  -'k. ' 


WVTFV7 


(34) 


J  -  £  /(£  +  -r  t1/3) 

max  o  o  2 


2  .  .  -1/3 

3  Co  +  T 


-  ,  -lr  1/3  -1  .  1/3,, 

i  t  3?  t  -a  arctan(o  t  )] 

1  4  O  o  O  j 

1  +  A  J  e  dr 

_ 2_ifl _ 

,  -l,~l/3  -1  .  ,  -1/3,, 

2-2/3,  35o  lT  “°o  arctan(aot  )] 
(1  r  a0T  )e 


The  generated  voltage,  UQ  *  RqI  +  L^dl/dt,  is  by  Eq.  (29)  in  dimensionless 
representation 

U(t)  -  J(t)  -  [(1  +  |  C0t"1/3)J(t)  -  |  50t"1/31M0(1  +  V2/3),  0  <  t  <  t. 

o 


The  inductive  singularity,  U( x)  -*■  for  t  -►  0,  is  due  to  the  divergence  of  the 
selfsimilar  speed  — ^ ,  d£(x)/dx  ■  (2Cq/3)t  +  »  for  r  -►  0,  and  has,  therefore, 

no  physical  meaning. 

-1  -1  5  -1 

Example.  For  b  *  10  m,  c  »  1  m,  x  *  10  m,  x  ■  1  a,  a  •  10  Q  /m, 

- * —  o  ® 

Rp  •  o-1c/b6  -  10_2n,  Ro  +  Rp  »  axuT1^  Lo  -  4irxl0-6  Henry,  v(10~6sec)  - 

4  -5 

2.650x10  m/sec,  B  *  ITesla,  the  generator  characteristics  are:  t  *  2irxl0  sec, 
o  o 

t  *  1.262xl0_4sec,  C  *  2ir/10,  A  -  11/10,  I  *  106/4ir  Amps,  I  »  lxlO4  Amps, 

o  o  o  max 

U  *  2x10 3  Volt,  (01)  *  2x10 7  Watt,  AE  -  Ult  -  103  Joule, 

max  max 

1/3 

Fig.  3  exhibits  the  current  J(x)  versus  aQx  for  *  2tt/10  with  Aq  *  l/4ir, 
10/4tt,  and  100/4ir  as  parameter  based  on  Eq.  (29).  In  each  case,  J(x)  rises  rap¬ 
idly  to  a  maximum  value  and  then  decays  slowly  with  increasing  time.  The  peak 

current  J  decreases  with  increasing  A  (inductance) *  It  should  be  noted  that 
max  o 

*1/3  1/2 

the  distributions  J(x)  are  applicable  up  to  the  transit  value  a^x  *  (x^/A^c) 

[Eq.  (28)]  which  depends  on  the  respective  system  parameters.  Experimental 

3) 

current  oscillograms  for  the  detonation  flow—  agree  only  qualitatively  with 
Eq.  (29)  since  the  electrodes  were  connected  in  their  centers  by  the  external 


load  circuit. 
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Current  transient  J(t)  of  detonation  shock  flow  versus 
1/2  1/3 

(C0/Aq)  t  for  £  ■  2tt/10,  and  A  ■  l/4ir;  10/4ir;  100/4* 


CONCLUSION 


An  elementary,  one-dimensional  theory  for  an  explosion  driven  in¬ 
duction  generator  has  been  developed,  in  which  the  upstream  ends  of  the 

electrodes  are  connected  by  an  external  R  -L  load  circuit.  The  formu- 

o  o 

las  derived  for  the  current  I(t)  and  voltage  U(t)  permit  an  optimum  gen¬ 
erator  design  by  appropriate  selection  of  the  system  geometry,  flow,  and 
load  parameters. 

The  jet  flow  with  (about)  constant  shock  speed  v(t)  *  vq  produces 

an  order-of-magnitude  more  energy  than  a  comparable  detonation  flow  with 

-1/3 

shock  speed  v(t)  «  t  decreasing  with  time.  Since  shock  speeds  in¬ 
creasing  with  time  are  not  realizable  in  the  coplanar  duct,  the  jet  flow 
is  ideal  for  power  generation. 

The  magnetic  self-field  B(t)  of  the  generated  current  I(t)  is  of 
the  same  order-of-magnitude  as  the  external  magnetic  field  Bq.  B(t) 
considerably  reduces  the  generator  performance  if  the  electrodes  are 
connected  at  their  upstream  ends  by  the  external  load  circuit  (l(t)  and 
1*o  antiparallel) . 

An  optimum  power  output  is  achieved  by  attaching  the  external  load 

circuit  to  the  downstream  ends  of  the  electrodes  so  that  B(t)  and  B  are 

o 

in  the  same  direction.  For  a  maximum  energy  output  the  parameters 
I q9  i  should  be  chosen  as  large  and  A^as  small  as  possible. 

For  the  transit  period  t  of  the  shockwave,  explosion  driven  gener¬ 
ators  are  generally  contained  by  inertial  confinement.  For  this  reason, 
no  attempt  has  been  made  to  discuss  generator  performance  for  times 
beyond  t. 
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DETONATION-DRIVEN  INDUCTION  GENERATOR 
WITH  POSITIVE  SUPERPOSITION  OF  EXTERNAL  AND  INDUCED  MAGNETIC  FIELDS 


H.  E.  Wilhelm 

Michelson  Laboratory,  Physics  Division 
Naval  Weapons  Center,  China  Lake,  California  93555 

ABSTRACT 

Based  on  Maxwell's  equations,  the  electric  current  and  voltage  pulses 

induced  in  explosion-driven  induction  generators  (R  «  u  ovL  »  1)  with  plane 

mo 

electrodes  and  homogeneous  external  magnetic  field  B^  are  calculated  in  closed 

form  for  an  external  load  circuit  with  resistance  R  and  inductance  L  ccm- 

o  o 

connected  to  the  downstream  ends  of  the  electrodes  so  that  the  magnetic  self- 
field  B(t)  of  the  generated  current  I(t)  is  in  the  same  direction  as  the 

-v 

transverse  external  field  Bq.  Two  plasma  shock  flow  models  are  considered: 

(i)  the  jet  flow  with  shock  speed  v(t)  *  vQH(t)  and  (ii)  the  plane  detonation 

flow  with  shock  speed  v(t)  «  (2/3)  (E  /p  )^^t  It  is  shown  that  the  jet 

o  o 

flow  produces  large  cumulative  current  and  voltage  pulses  which  increase  with 
time,  whereas  the  detonation  flow  generates  considerably  smaller  and  shorter 
current  and  voltage  pulses.  In  both  cases,  the  magnetic  self-field  B(t)  of  the 
generated  current  I(t)  is  of  the  order  of  magnitude  of  Bq.  In  case  of  positive 
superposition  of  the  Induced  magnetic  field  on  the  external  magnetic  field  Bq 
maximum  flux  changes  are  obtained.  For  this  reason,  in  induction  generators 
the  external  load  circuit  should  not  be  attached  to  the  upstream  ends  or  the 
centers,  but  to  the  downstream  ends  of  the  electrodes. 
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INTRODUCTION 


Induction  generators  are  based  on  the  electromagnetic  interaction  of  an 
electrically  conducting  flow  field  with  a  magnetic  field,  where  at  least  one  of 
the  fields  is  transient  so  that  voltages  are  induced  in  an  external  circuit 
which  are  due  to  temporal  magnetic  flux  changes  d$/dt  ^  0.  An  explosion-driven 
induction  generator  with  a  transient  plasma  shock  flow  in  a  constant  external 
magnetic  field  Bq  has  been  analyzed^  when  the  external  Rq  -  Lq  load  circuit  is 
attached  to  the  flow  entrance  ends  of  the  electrodes.  In  this  case,  the  self- 
magnetic  field  B(t)  of  the  generated  current  I(t)  is  antiparallel  to  the  trans¬ 
verse  external  magnetic  field  Bq  so  that  the  magnetic  flux  changes  d$/dt  during 
the  transit  time  t  of  the  flow  are  reduced.—^  Herein,  a  novel  explosion-driven 
induction  generator  with  optimum  magnetic  flux  changes  is  proposed,  in  which 

the  external  R  -  L  load  circuit  is  connected  to  the  flow  exit  ends  of  the 
o  o 

coplanar  electrodes  so  that  the  induced  magnetic  field  B(t)  of  the  current  I(t) 

in  the  shock  front  current  layer  and  the  electrodes  can  superimpose  themselves 

-*• 

positively  on  the  external  magnetic  field  Bq.  This  effect  is  quantitatively 
significant  since  B(t)  and  Bq  are  parallel  and  of  the  same  order  of  magnitude. 
Thus,  considerably  larger  current,  voltage,  and  power  pulses  are  induced  than 
in  the  antiparallel  Bq  -  B  case. 

In  induction  generator  experiments,  the  external  Rq  -  Lq  load  circuit  is 

2  3  4) 

usually  attached  to  the  center  of  the  electrodes.  *  f  In  earlier  experiments, 

a  detonation-produced  plasma  with  a  decreasing  shock  velocity  v(t)  0  for  t  -►  • 

2) 

is  directly  injected  into  the  inter-electrode  space,—  whereas,  in  more  recent 
3.4) 

experiments,—1—  discontinuous  plasma  jets  are  driven  across  an  external, 
constant  magnetic  field  Bq.  Plasma  jets  v(t)  *  vQH(t)  of  nearly  constant  shock 
speed  vq  are  produced,  e.g«,  by  implosion  of  a  dense  noble  gas  between  flyer 
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3 ) 

plates  and  subsequent  expansion  through  a  diaphragm  into  the  generator  channel,— 

or  by  means  of  a  tube-like  arrangement  of  the  solid  explosive  resulting  in  a 

cylindrical  implosion  with  large  axial  shock  velocity  of  the  filling  gas.— ^ 

Detonation-driven  induction  generators  have  been  operated  with  shock  veloci- 
A  5 

ties  v  *  10  -  10  m/sec,  external  magnetic  fields  B  »  1  -  10  Tesla,  conductivi- 

A  “1  1  2  10  12  —3 

ties  a  >  10  0  /m,  and  flow  energies  y  pv  "10  -  10  Joule  m  at  pressures 

A  2-3 

p  >  10  bar  (p  >  10  kg  m  ) .  The  corresponding  magnetic  Reynolds  number  R 

m 

is  large,  whereas  the  magnetic  interaction  number  R^  is  small: 

R  =  u  ovL  »  1  ,  R .  -  B2/u  pv2  «  1  ,  (1) 

mo  loo 

where  L  -  |Vxb/b|.  For  Rffl  »  1,  the  electric  current  is  restricted  to  a  thin 

layer  6  ~  (t/uQa)  behind  the  shock  front  (t  »  transit  time  of  flow  in  Bq).  In 

1  2 

view  of  the  large  dynamic  pressure  y  pv  compared  to  the  magnetic  pressure 
2  6  8  —3 

B*/2pq  *  10  -  10  Joule  m  ,  the  efficiency  of  energy  transformation  is  small, 

while  the  absolute  performance  is  considerable  since  currents  I  -  (B  /u  )b  > 

o  o 

A  —a  A 

10  amps  (b  *  duct  width)  and  energies  Ult  >  10  Joule  have  been  obtained  experi- 

3  A) 

mentally  per  detonation.—2— 

In  the  following,  the  initial- value  problem  for  a  detonation -driven  generator 

is  solved  in  closed  form  by  means  of  Maxwell’s  equations,  in  which  the  downstream 

ends  (x  *  x  )  of  the  coplanar  electrodes  are  connected  by  an  external  L  -  R 
o  r  y  o  o 

load  circuit.  In  this  case,  solutions  for  the  magnetic  field  B(t)  of  the 
generated  currents  I(t)  exist  which  are  parallel  to  the  external  magnetic  field 
Bq  (Fig.  1).  Two  types  of  plasma  shock  flows  are  considered  corresponding  to 
i)  the  idealized  jet  flow  with  constant  shock  speed  vq  and  il)  the  plane  detona¬ 
tion  flow  due  to  an  explosive  energy  release  Eq  per  unit  area  in  a  gas  of 
initial  density  p^: 

v(t)  -  voH(t)  ,  x(t)  -  vQtH(t)  ,  (2) 

0  <  t  <  t  ,  t  ■  x,/v0  » 
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field  i  ,  and  R  -  L  load  circuit  CD.  Shock  front  AB  at 

r>  '  a  A 


and 


v(t)  =  (2/3)(Eo/po)1/3t"1/3  ,  x(t)  =  (Eo/po)1/3t2/3  ,  (3) 

0  <  t  <  t  ,  t  =  (Eo/po)-1/2xo/2 

where  x(t)  is  the  shock  front  position,  and  t  is  the  transit  time  of  the  shock 
front  x(t)  in  the  duct  of  effective  length  xq.  The  self -similar  solution  (3) 
is  asymptotically  correct  for  larger  times  t ^  but  can  be  used  for  t  >  0  since 
the  singularity  of  v(t)  at  t  -  0  is  without  physical  consequences  for  generator 
applications. 

The  mechanical  stability  of  the  generator  during  the  shock  transit  time  t 
is  provided  for  by  inertial  confinement  of  the  duct  through  high  density  plaster. 
No  attempt  is  made  at  discussing  the  electrical  output  for  times  t  >  t  of  the 
explosive  disintegration  of  the  generator.  The  external  circuit  location  xq  is 
an  effective  one,  which  can  be  calculated  once  the  geometries  of  the  inductance 
L  and  resistance  R  are  specified. 
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INITIAL -VALUE  PROBLEM 


Figure  1  exhibits  a  one-dimensional  model  for  a  detonation-driven  generator 

in  which  the  shock  wave  flow  v(t)  =  dx(t)/dt  is  initiated  at  time  t  s  0  in  the 

plane  x  =  0  by  explosive  charges  in  the  space  x  <  0  so  that  the  shock  front  is 

at  x  =  x(t)  >  0  for  t  >  0.  The  copper  electrodes  in  the  planes  z  *  ±c/2  have 

the  extensions  x  >  xq  <  00  and  -b/2  <  y  <  +b/2  and  are  connected  by  an  external 

circuit  CD  with  load  R  and  inductance  L  at  the  flow  exit  end  x  =  x  .  The 

o  o  o 

entire  system  is  embedded  in  a  homogeneous  magnetic  field  Bq  *  {0,Bq,0}.  As  a 

result  of  the  v  x  Bq  interaction  of  the  conducting  shock  front  (a  >  10^  ft  ^/m) 

with  the  external  magnetic  field  Bq,  a  current  I(t)  is  induced  behind  the  shock 

*  1/2 

front  AB  in  a  cross  section  AxAy  =*  b6,  where  a  ~  (t/uQa)  is  the  (small)  width 

of  the  current  layer.  The  current  I(t)  flows  in  the  closed  loop  A  B  C  ■+ 

D  A,  where  the  leg  AB  at  x  *  x(t)  moves  with  the  shock  speed  v(t)  *  dx(t)/dt. 

Application  of  the  surface  curl  equation  n  x  [B]  =  y  I(t)/b  to  the  inter - 

o 

faces  AB,  BC,  CD,  DA  demonstrates  —  in  the  one-dimensional  approximation 
neglecting  end  effects  at  y  =  ±b/2  and  x  *  xq  —  that  the  current  I(t)  produce? 
a  magnetic  self-field  in  the  extending  loop  volume, 

B(t)  *  {0,woI(t)/b,0}  ,  x(t)  <  x  <  xq  ,  |y|  <  y  ,  |z|  <  -|  ,  (4) 

which  is  parallel  to  the  external  magnetic  field  if  I(t)  >  0.  The  total 
magnetic  field  within  the  loop  volume  is  the  positive  superposition 

B(t)  -  {0,Bo+uoI(t)/b,0}  ,  x(t)  <  x  <  xQ  ,  |y|  <  j  ,  |z|  <  j  ,  (5) 

whereas  $(t)  ■  Bq  outside  of  this  region.  Integration  of  the  Maxwell  equation 
v  x  1  -  -a!/3t  over  the  cross  section  ABCD  of  the  extending  current  loop  yields 
by  Stokes'  theorem 
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dr  *  -  $  *  dS  +  <p  (v  x  b)  *dr 


(6) 


where 


E  *  E  +  v  x  B 


«  1/p  e 

o  o 


(7) 


and  v(t)  is  the  velocity  field  of  the  sections  of  the  closed  loop,  at  the  point 
corresponding  to  the  path  element  dr  [v(t)  =  exdx(t)/dt  for  AB,  but  v(t)  *  0  for 
CD,  DA,  CB] .  The  surface  element  d!s  =  +e^cdx  forms  with  dr  [pointing  everywhere 
in  the  direction  a  right-hand  system.  By  Eqs.  (6)  and  (7) 

«o  +  VI(t)  +  Lo  ■  ->=  C°  d*  +  I*'2  +  *<‘>1*  •  <«> 

x(t)  -c/2 

B  *  D 

since  /  E  dz  31  IRp  (Rp  »  plasma  resistance),  J  Edz  *  RqI  4*  L^dl/dt,  and  / Edz  *  0 

^  C  g  —  ^ 

for  the  sections  BC  and  DA  of  the  highly  conducting  electrodes  (a  *  10  /m) . 

v^u 


Accordingly, 


<*„  +  VI(t)  +  3t  {tlo  +  £«=)IICt:)}  -  cBo  ^ 


(9) 


where 


L(t)  *  pqc[xo  ~  x(t)]/b  >  0  , 


0  <  x(t)  <  x 


(10) 


is  the  time-dependent  inductance  of  the  electrode  sections  AD  and  BC  of  decreasing 
length  Ax  -  xq  -  x(t).  Multiplication  of  Eq.  (9)  by  I(t)  leads  to  the  energy 
conservation  law 


d  rl 


(Ro  +  RP)I  +  3t  {I  tLo  +  L(t)]l‘} 


UoC  t2  dx(t) 


2b 


dt 


♦'TT1*.' 


(11) 


The  plasma  resistance  Rp  *  Rp(t)  is  a  slowly  varying  function  of  t  and  can 
be  calculated  only  via  a  detailed  physical  model  of  the  current  sheath  5  at 
x  *  x(t).  However,  since  Rq  is  sufficiently  large  in  actual  experiments, 

Rq  +  Rp  can  be  treated  as  a  quasi-constant. 
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Rq  +  Rp(t)  =  constant  ,  Rq  >  Rp  ;  Rq  +  Rp(t)  =  Rq  ,  Rq  >>  Rp  •  (12) 

For  a  generalized  formulation,  dimensionless  dependent  and  independent 
variables  are  introduced  by 

J(t)  =  I(t)/IQ  ,  T  =  t/tQ  ,  c(x)  =  x(t) / c  (13) 

with 

Io  *  (Bo/uo)b  ’  Co  "  <V2/b>/(R0  +  V  »  Ao  =  (Vc)  +  Lo^V2/b)  *  (14) 

By  Eqs.  (9),  (10),  (13),  and  (14),  the  electric  current  I(t)  in  the  explosion- 
driven  generator  with  positive  external  and  induced  magnetic  field  superposition 
is  determined  by  the  inhomogeneous  initial-value  problem  with  variable  coeffi¬ 
cients: 

£  {[AQ  -  l(T)]J(T)}  +  J(t)  -  ,  0<T<T  ,  (15) 

J(t  ■  0)  =  0  .  (16) 

Equations  (15)  and  (16)  are  solved  for  the  jet  and  detonation  flows  given  by 
Eqs*  (2)  and  (3),  respectively,  which  define  the  dimensionless  coordinate  £(t) 
of  the  shock  front. 

12  2 

In  view  of  the  large  flow  energy,  j  pv  >>  Bq/2uo,  velocity  decreases  due 
to  Lorentz  and  viscous  forces  are  neglected  in  the  above  generator  model.  Further¬ 
more,  viscous  boundary  layers  are  disregarded  since  the  plasma  moves  in  the  form 
of  a  rectangular  velocity  slug  through  the  duct  In  large  Reynolds  number  flow, 

R  ■  pvL/y  >>  1. 
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For  the  jet  flow,  the  velocity  and  coordinate  of  the  shock  front  are  given 

in  Eq.  (2).  The  dimensionless  shock  front  coordinate  €(t),  £  ,  and  the  transit 

o 

time  t  are,  by  Eqs. (13)  and  (2), 


,  +Rp 

“5— * 1  )  •  (17) 

c  /b 


V  u  c  /  b 

-  t  (  rVr;  > 


x  R 

-S.  (  - 

v 

O  u 


(23) 


Jmax  1  5o/(1  "  V  for  <  1  ,  Lq  «  UQcxo/b  , 


Jmax  "  -  *n(1  -  —  x)  <  «  for  ?o  «  1  . 

o 


(24) 


The  generated  voltage,  U  *  RqI  +  L  dl/dt*  is  by  Eq.  (18)  in  dimensionless  repre¬ 
sentation  (U  =  i  R  ) 
o  o  o 

<%>  -r^r  [1  -  o  -  rT>a'5»)/s°i  o  -rT)<1'25°)/E»  • 

o  o  o  o  o 


0  <  T  <  x  ,  (25) 


where 


11m  U( t)  *  -  £n(l  -  t) 


L  /R 


o_  o  .-1  .-1  ,-l 


Co  -  1 


t  Aq  (1  -  A0  t)  "■<-,  0<T<T 


(26) 


and 


U(t)  -  [1  +-t— 5 "  ( 


L  /R  1  -  2£  5  L  /R  C  § 

oo,  °  m  °  r  i  °  °  °  o 

A -  )]A“T  ~  T  *  A-  T  1  * 

o  o  o  o  o 


(27) 


U 


-  t-V  u  -  a  -  ^  a  -  r  f)(1-2!<.)/5<.  . 

max  1  -  £  1  v  A  7  t  A  A  7  * 

o  o  o  o  o 

T  =  T  .  (28) 


For  special  cases,  the  maximum  voltage  is  given  by  simpler  formulas, 


max 


V(1  “9  for  5o<  1  ’  Lo  ^  uocxo/b  * 


,  L  /R  1  -1  -1 

U  m  -  £n(l  -  A  t)  +  —7 — —  A  ^(l  -  A  x)  <  ®  for  ”  1 
max  v  o  t  o  v  o'  o 


(29) 

(30) 


Example.  For  b  ■  10  ^  m,  c  *  1  m,  xq  *  1  m,  o  ■  10^  0  Vm,  Rp  a  o  ^c/b6  ■ 

10*2  0,  Rq  +  Rp  -  2  x  10_1  0,  Lq  -  4ir  x  10*6  Henry,  vq  -  104  m/sec,  Bq  - 

1  Tesla,  the  characteristics  of  the  generator  are:  tQ  ■  2ir  x  10  ^  sec,  t  ■ 

10"4  sec,  -  2ir/10,  Aq  -  2,  lQ  -  106/4ir  amps,  1^  s  1  x  105  amps,  Umt  a 

2  x  104  volts,  (UI)  a  2  x  109  watts,  AE  -  Olt  -  105  joules. 

*  max  * 
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Figure  2  shows  the  current  J(x)  versus  (So/Aq)t  for  ■  2ir/10  and  * 

9/10  in  the  case  of  L  <<  y  cx  / b  so  that  t  g  t  based  on  Eq.  (18).  While  for 

o  o  o 

£  =  2tt/10  the  current  rises  steadily  to  J  *  2tt/(10  -  2ir)  *  1.690,  J(x)  *  1 

o  max 

nearly  throughout  the  interval  0  <  t  <  x  for  *  9/10,  and  the  sharp  rise  to 

J  *  9  occurs  not  until  the  very  end  t  £  x  of  the  flow  transit  period.  It 
max 

should  be  noted  that  the  generator  produces  in  each  case  average  currents 

T  >  I  in  excess  of  the  normalization  current  X  *  (B  /y  )b.  Experimental 
o  o  o  o 

current  oscillograms  are  not  available  for  comparison,  since  the  advantages  of 
a  generator  with  a  load  circuit  attached  to  the  downstream  ends  of  the  elec¬ 
trodes  have  apparently  not  yet  been  recognized. 
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J 

max 


*  «0'tf<w1/2  -  <„> f"  «„  *  (i)2/3»y3 


L  «  m  cx  /b 
o  o  o 


(37) 


Eq.  (32)  has  positive,  monotonically  increasing  solutions  J(x)  >  0  for  t  c(0,t). 
The  magnetic  field  l(t)  of  this  current  is  in  the  same  direction  as  the  external 
magnetic  field  Bq  (positive  superposition) . 

The  generated  voltage,  U  *  RqI  +  LQdI/dt,  is  by  Eq.  (32)  in  dimensionless 
representation  (Uq  -  IqRo)  * 


L  /R 

U(t)  *  J(t)  +  ~  —  f 


(2/3)£0  -  (t 


1/3 


-  (2/3U  ]J(t) 
o  1 


A  2  2/3N  1/3 

Ao(1  •  v  )T 


0  <  T  <  T 


(38) 


The  inductive  singularity r  U(t)  00  for  r  0,  is  due  to  the  divergence  of  the 
self-similar  speed,— ^  JE(r)/dT  -  (2£q/3)t  +  m  for  t  0,  and  has,  therefore, 

no  physical  meaning. 

Example.  For  b  »  10  ^  m,  c  *  1  m,  x  »  1  a,  0  »  10^  ft  ^/m,  R_  *  a  ^“c/b<5  * 
- - -  o  p 

10~2  Q,  Ro  +  Rp  *  2  X  lO-1  fl,  Lo  »  4ir  X  lO-6  Henry,  v(t  -  10~6  sec)  -  2.650  x 

4  -5 

10  m/sec,  Bq  ■  1  Tesla,  the  generator  characteristics  are:  tQ  *  2tt  x  10  sec, 

t  *  1.262  x  10  4  sec,  £  *  2tt/10,  A  *  2,  I  *  10^/4*  amps,  I  *  3xl04  amps, 

o  o  o  max 

U  *  6xl03  volts,  (UI)  *  2xl08  watts,  AE  -  Olt  -  104  joules, 
max  *  max  J 

1/2  1/3 

Figure  3  shows  the  current  J(t>  versus  (C  /A  )  t  for  £  *  2ir/10, 

o  o  o 

-  9/10,  and  Aq  *  2  based  on  Eq.  (32).  In  each  case,  J(t)  increases  monoto- 

1/2  1/3 

nically  from  0  to  its  final  value  J  at  transit  time  t,  where  (5  /A  )  '  x 1  * 

max  o  o 

1/2 

(x  /A  c;  •  Jl/2  for  x  »  1  m,  c  *  1  m,  and  A  -  2.  On  the  average,  the  pro- 
o  o  o  o 

duced  currents  f  <  I  are  smaller  than  the  normalization  current  I  ■  (B  / u  )b. 

o  o  o  o 

Comparison  with  Fig.  2  indicates  that  J(t)  of  the  detonation  flow  is  by  an  order 
of  magnitude  smaller  than  J(t)  of  the  jet  flow. 
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CONCLUSIONS 


An  electromagnetic  theory  for  an  explosion-driven  induction  generator  has 

been  developed  in  which  the  downstream  ends  of  the  electrodes  are  connected  by 

an  external  R  -  L  load  circuit.  The  formulas  derived  for  the  current  I(t) 
o  o 

and  voltage  U(t)  permit  an  optimum  generator  design  by  appropriate  selection 
of  the  system  geometry,  flow,  and  load  parameters. 

The  jet  flow  with  (about)  constant  shock  speed  v(t)  £  vq  produces  consid¬ 
erably  more  energy  than  a  comparable  detonation  flow  with  shock  speed  v(t)  *  t~^ 
decreasing  with  time.  Since  shock  speeds  increasing  with  time  are  not  realizable 
in  the  coplanar  duct,  the  jet  flow  is  ideal  for  power  generation. 

The  magnetic  self-field  B(t)  of  the  generated  current  I(t)  is  of  the  same 
order  of  magnitude  as  the  external  magnetic  field  B^.  The  induced  magnetic 
field  ?(t)  considerably  increases  the  performance  of  generators  with  shock 
flow,  since  it  superimposes  itself  positively  on  the  external  magnetic  field  Bq 
if  the  electrodes  are  connected  at  their  downstream  ends  by  the  external  load 
circuit  (B(t)  and  Bq  in  the  same  direction  for  optimum  magnetic  flux  changes). 

An  optimum  power  output  is  achieved  by  (i)  attaching  the  external  load 
circuit  to  the  downstream  ends  of  the  electrodes  and  (ii)  driving  the  generator 
with  a  jet  shock  flow.  For  a  maximum  energy  output,  the  parameters  Iq,  t 

should  be  chosen  as  large  as  possible,  and  should  be  chosen  as  small  as 
possible. 
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POWER  GENERATION  BY  SHOCKWAVE  INTERACTION  WITH  SOLIDS 


ELECTRICALLY  POLARIZABLE  BY  STRESS 

H.E.  Wilhelm 

Michelson  Laboratory,  Naval  Weapons  Center,  China  Lake,  CA  93555 


Abstract 

An  elementary  theory  of  the  interaction  of  shock  waves  with 
solids  which  are  electrically  polarizable  by  stress  is  given  and 
discussed  with  regard  to  the  generation  of  electric  power  pulses. 
The  solid  is  assumed  to  be  a  one-dimensional  slab  between  plane 
electrodes,  which  are  connected  through  an  external  load  circuit. 

It  is  shown  that  the  temporal  current  build-up  in  the  external  cir¬ 
cuit,  due  to  electric  polarization  and  electric  conduction  behind 
the  shock  front  of  the  stress  wave,  space  charge  and  electric  po¬ 
larization  relaxations,  is  determined  by  an  inhomogeneous  Volterra 
integral  equation.  The  latter  is  solved  analytically,  and  current 
pulse  forms  for  incident  step-shaped  stress  waves  are  calculated. 
For  shocked  solids  with  cross  section  of  the  order  of  one  square 
meter,  energies  of  the  order  AE  -  10^  Joule  per  pulse  should  be 
achievable. 


37 


INTRODUCTION 


Although  electric  energy  pulses  of  10  Joules  and  larger  have  been 

1-3) 

produced  in  explosion  driven  induction  generators  -  ,  the  requirement 

of  large  external  magnetic  fields  of  at  least  1  Tesla  renders  these  mag¬ 
neto  gas  dynamic  generators  unsuitable  for  many  applications,  in  particular 
since  the  prospects  for  self -excitation  of  such  magnetic  fields  during 
pulse  times  of  the  order  10  ^  3ec  are  not  promising  .  For  these  reasons, 
other  methods  for  conversion  of  chemical  energy  of  explosives  into  elec¬ 
tric  energy  have  to  be  found.  E.g.,  such  a  "non-magnetic"  power  supply 
could  be  used  to  produce  the  magnetic  field  required  for  magnetogas dyna¬ 
mic  generators. 

Not  only  piezoelectric  crystals  such  as  tourmaline  or  ferroelectric 

crystals  such  as  lead  circonate  titanate,  but  also  a  large  number  of 

other  solids  become  strongly  electrically  polarized  when  subject  to  pres¬ 
et 

sures  p  >  10  bar.  Such  materials  are,  e.g.,  electrets,  semiconductors, 

4) 

alkali  halides,  plexiglas ,  polysterene  —  .  These  solids  simultaneously 

develop  good  to  metallic  electrical  conductivities  when  subject  to  pres- 
4  6 

sure  p  *  10  -10  bar,  due  to  non thermal  ionization  by  overlapping  of  at¬ 
omic  wave  functions  at  extreme  pressures  — ^ .  Similar  to  the  chemical 
polarization  field  in  a  battery,  the  electric  polarization  caused  by  the 
stress  shock  wave  in  the  solid  drives  an  electric  current  pulse  through 
an  external  circuit  attached  via  electrodes  (transverse  to  the  polariza¬ 
tion  vector)  to  the  solid. 

The  interaction  of  a  shock  wave  with  solids,  which  are  electrically 

polarizable  by  stress,  appears  to  have  potential  for  the  generation  of 

electric  power  pulses.  Table  I  gives  the  stress  polarization  coefficient  — ^ 

f,  the  dielectric  permittivity  ^  e  -  e*eo>  the  density  p,  and  (dimensional 
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TABLE  I:  Order-of-magnitude  estimates  of  stress  polarization  P,  electric 

2 

field  E,  current  I,  and  energy  AE  for  a  stress  polarizable  slab  (A  *  1  m  , 

-2  3 

a  *  10  m)  subject  to  a  shock  wave  with  speed  v  *  10  m/sec. 


Quartz 

BaTi03 

Rochelle  Salt 

f [A  sec/N] 

-12 
2.3x10  J' 

l.OxlO-10 

A.OxlO-10 

p[kg  m  3] 

2.65xl03 

6.02xl03 

1. 77xl03 

o [N/m2] 

1.3xl09 

3. lxlO9 

0. 9xl09 

«.  [A  sec/Vm] 

4x 10-11 

2xl0'8 

9xl0-11 

2 

P[A  sec/m  ] 

3.0xl0'3 

3. lxlO-1 

3.  lxlO-1 

E[V/m] 

7.6xl07 

1.5xl07 

3.9xl09 

I  [A] 

3xl02 

3xl04 

3xl04 

AE[J] 

lxlO5 

5xl06 

lxlO8 
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*12 

analysis)  order- of -magnitude  estimates  of  the  stress  amplitude  a  *  ypv  , 

electric  polarization  P  -  af,  electric  field  E  -  P/c  ,  short  circuit  current 

2 

I  -  APv/a*  and  electric  energy  AE  -  (P  /2e)Aa  for  one  piezoelectric  and  two 

2  -2 

ferroelectric  slabs  of  cross  section  A  *  1  m  and  thickness  a  *  10  m  in  case 

3 

of  a  stress  wove  with  (moderate)  speed  v  *  10  m/sec.  It  is  seen  that  E 

of  quartz  is  of  the  order  of,  E  of  BaTiO^  is  below,  and  E  of  Rochelle  salt 

4  ) 

is  above,  the  dielectric  breakdown  field  —  .  While  the  generated  voltages 

2 

V  -  Ea  are  large,  the  generated  currents  I  over  an  area  A  ■  1/m  are  only 
moderate,  but  the  energy  pulses  AE  are  significant.  Furthermore,  these 
examples  show  that  a  stress  polarization  generator  can  be  operated  both 
below  and  in  the  breakdown  regime. 

The  electric  polarization  of  solids  by  stress  shock  waves  is  also  of 

considerable  basic  physical  interest.  Attempts  at  a  shock  wave  interaction 

theory  with  stress  polarizable  solids  (short  circuited)  were  made  by  Alli- 
7)  8)  9) 

son  —  ,  Zeldovich  — '  >  and  Ivanov  —  under  various  phenomenological  assump¬ 
tions  and  compared  with  experimental  data  from  piezoelectric  and  ferroelec- 

4) 

trie  pressure  transducers  —  .  The  applicability  and  limitations  of  these 
theories  have  been  the  subject  of  extensive  discussions  ,  which  in¬ 

dicate  that  the  problem  is  not  completely  understood  from  first  principles. 

Herein,  a  theory  of  the  interaction  of  stress  shock  waves  with  (non- 
ferroelectric)  stress  polarizable  solids  is  presented  based  on  Maxwell's 
equations,  when  the  external  circuit  contains  an  ohmic  load.  It  is  shown 
that  the  electric  current  pulse  I(t)  is  determined  by  an  inhomogeneous  Vol- 
terra  Integral  equation  of  convolution  type.  Closed-form  solutions  are  ob¬ 
tained  for  the  generated  current  I(t)  by  taking  into  consideration  the  dif¬ 
ferent  time  scales  at  which  electric  space  charge  and  electric  polarization 
by  stress  develop  in  the  solid. 
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The  following  shock  wave  interaction  theory  with  solids  is  in  general 
not  applicable  to  ferroelectric  solids  (except  for  rough  estimates).  Ferro- 
electrics  exhibit  large  permanent  electric  polarization  and  a  pronounced 

"V 

hysteresis  effect  so  that  the  electric  displacement  D  *  D(E)  is  no  longer  a 
linear  but  a  strongly  nonlinear  function  of  the  electric  field  E  (a  hyster¬ 
esis  theory  for  ferroelectrics  is  not  available) . 
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POLARIZATION  AND  CHARGE  BUILD-UP 

In  unstressed,  shock  polarizable  (non-ferroelectric)  solids,  the  electric 
dipole  moments  of  the  molecules  add  up  to  a  zero  vector,  whereas  under  stress 
the  dipoles  are  displaced  preferentially  in  one  direction  so  that  a  net  polar¬ 
ization  results*  This  is  illustrated  in  Fig.  la  (a  *  0)  and  Fig.  lb  (a  4  0) , 

which  depicts  the  unstressed  and  stressed  dipole  configurations  in  a  quartz 
crystal  subject  to  a  stress  component  o  which  is  perpendicular  to  the  optical 
axis.  In  Fig.  lb,  points  to  the  left,  but  ?  would  point  to  the  right  if 
the  crystal  were  rotated  by  180°.  For  this  reason,  the  polarization  vector 
*£  may  point  in  or  opposite  to  the  direction  of  propagation  of  the  stress  wave 
depending  on  the  chosen  orientation  of  the  crystal.  For  a  stress  wave  o(x,t) 
propagating  in  the  positive  x-direction,  the  designation  ,f±"  polarization  is 
used  if  points  into  the  +x-directions.  The  resulting  direction  of  the  cur¬ 
rent  I(t)  in  the  external  circuit  is  shown  in  Fig.  2a  for  the  "+"  polariza¬ 

tion,  and  in  Fig.  2b  for  the  polarization,  based  on  Kirchhoff’s  law 
j  l*dr  -  0. 

The  stress  polarization  generator  model  (Fig.  3)  consists  of  a  shock 
polarizable  solid  of  cross  section  A  and  thickness  a  between  plane  electrodes 
at  x  ■  0  and  x  *  a,  which  are  connected  through  an  external  circuit  with  load 
R  >  0.  At  time  t  -  0,  an  explosion  produced,  step  shaped  stress  wave  o(x,t) 
with  amplitude  a  and  propagation  speed  v  impinges  normal  onto  the  surface 
x  ■  0  of  the  solid  so  that,  at  time  t  >  0,  the  shock  front  is  at  x  ■  vt  and 
the  stress  distribution  in  0  <  x  <  a  is 

o(x,t)  -  oH(vt  -  x)  ,  t  >  0  ,  (1) 
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where 


H(vt  -  x)  *  1 
0 


0  <  x  <  vt 
vt  <  x  <  a 


,  (2) 


7  8  9) 

is  the  Heaviside  function  (compression  of  the  slab  Ax  ■  a  neglected)  7  2  J 
Behind  the  shock  front  x  -  vt,  the  polarization  P(x,t)  builds  up  according 
to  a  relaxation  process  with  characteristic  time  x, 


with 


P(x,t  »  0)  *  0 


(3) 

(4) 


as  initial  condition  for  the  ,,±”  polarizations  (Fig.  2a, b),  since  the  solid  is 
initially  unpolarized  and  unstressed.  In  view  of  Eq.  (1),  the  solution  of 
Eqs.  (3)  -  (4)  is 


P(x,t)  -  ±fa [1  -  e'<vt'x)/vT]H(vt  -  x)  .  (5) 

Accordingly,  P(x,t)  drops  exponentially  to  zero  for  x  +  vt  (shock  front),  and 
P(x,t)  *  ±fa  in0<x<vt  for  t  >>  t.  The  polarization  relaxation 
time  is  of  the  order-of-magnitude  t  -  10  sec.— 

Due  to  shock  polarization  P(x,t)  in  the  stressed  region,  which  is  absent 
in  the  unstressed  region,  the  electric  displacements  behind  and  ahead  of 
the  shock  front  x  *  vt  are  — ^ 


elEl  +  p  * 

0  <  x  <  vt 

> 

(6) 

*2®; 2  * 

vt  <  x  <  a 

9 

(7) 
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build-up  in  the  stressed  region  is  determined  by 


_of_  e-(vt-x)/vrH(vt  _  x) 

T  VT 
O 


0  <  x  <  vt 


(14) 


Px(x,t  »  0)  -  0 


(15) 


Hence 


Px(x,t) 


-  qf 

+  v(T  -  TQ) 


r  -(Vt-x) /VT 

[e 


-(vt-x)/vT 

-e  ]H(vt  -  x) 


(16) 


E1(x,t)  -  E^(vt,t)  -  (vT{)/el)p1(x,t)  ,  (17) 

for  the  "±"  polarization  cases.  The  space  charge  relaxation  time  is  defined 
by 


to  -  e1/° 

Substitution  of  Eq.  (17)  into  Eq.  (9)  shows  that  the  x- integration 
E^(vt,t)  is  determined  by 

dE^(vt,t)  E^vt.t) 

dt  +  t  "  Ae 

o  1 

E^O.O)  -  0 


(18) 


"constant" 


,  (19) 

•  (20) 


Hence 

.  -t/r  t  +t'/x 

E.(vt,t)  -  (Ae.)-1*  7  I(t’)e  °dt’  .  (21) 

A  0 


Thus,  we  obtain  the  electric  fields  behind  and  ahead  of  the  shock  front  as  in¬ 
tegral  functionals  of  I(t): 


-t/r 


o  t  +t'/t 

-/ 
l  o 


E1  <1.0-^/  Kf>. 


of  ,  -<«-x>v, 


-  e 


(22) 


0  <  x  <  vt 
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vt  <  x  <  a 


(23) 


,  t 

E-(x,t)  *  (Ae  )  /  I(t')dt' 

1  0 

for  the  "±"  polarization  cases.  Since  =  0  in  the  unstressed  region, 

Eq.  (23)  obtains  by  direct  integration  of  Eq.  (10). 

The  surface  charge  density  p*(vt,t)  *  n*[D]x-vt  of  the  shock  front  is  by 
Eqs.  (5),  (6),  (7),  (22),  and  (23) 

i  c  1  "t/Tn  t  +t'/Tn 

p*(vt,t)  -  (Ae.)"1/  I(t’)dt’  -  (Ae  )_1e  °J  I(t’)e  °dt'  .  (24) 

L  0  1  0 

It  is  obvious  that  in  general  p*(vt,t)  +  0,  i.e.,  p*(vt,t)  ■  0  only  if  e^  *  e^ 
and  tq  -►  0.  Therefore,  the  assumption  D^(x,t)  ■  D2(x,t)  for  x  *  vt  made  in  the 
theory  of  piezoelectric  and  ferroelectric  pressure  transducers  is  in  general 
incorrect. 
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INTEGRAL  EQUATION  FOR  I(t) 


The  electric  current  pulse  I(t)  in  the  external  circuit  is  determined 
by  Kirchhoff 's  law  y  E-dr  *  0  for  the  closed  current  loop  with  load  R  in 
Fig.  3,  where  E(x,t)  in  the  solid  is  given  by  Eqr .  (22)  and  (23)  and  /£R*dr  * 
+  IR.  Thus,  we  obtain  for  the  transient  current  I(t)  the  inhomogeneous 
Volterra  integral  equation: 


vt 

Ae 


t  -(t-t’)/To 

J  e 

1  0 


I(t')dt’  + 


SfVT 

o 


t 

/  IftOdt’  +  RI(t) 
0 


-t/x  -t/To 
te  -re 

<*“ - TT1 - ■> 


(25) 


The  upper/lower  signs  refer  to  the  "±"  polarization  cases.  A  more  general  re¬ 
presentation  of  this  integral  equation  is  achieved  by  introducing  the  dimension¬ 
less  variables 

J(t)  -  I(t)/Io  ,  I0  =  a (cf/e^)A  ;  t  -  t/tQ  ,  tQ  =  t  5  (26) 

and 

d  *  a/vt  ,  R  »  R/(c  ^vr/A)  ,  e  *  e]/e2  *  6  *  TJX  *  (27) 

as  characteristic  dimensionless  parameters  (o  ^vx/A  *  internal  resistance  of 
solid  slab  of  thickness  vt  traversed  by  shock  during  polarization  relaxation 
time  x).  By  means  of  Eqs.  (26)  and  (27),  the  determination  of  J (£)  is  reduced 
to  the  initial-value  problem: 

f  +  Rj(t) 

0  0 

tf  *  -*£  /6 

-  ±6  (1  -  -  - - )  ,  (28) 
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■*&*$?’* 


J(£«0)  *  0 


(29) 


By  means  of  the  Laplace  transformation  and  the  convolution  theorem,  the  integral 
equation  (28)  is  reducible  to  an  inhomogeneous  differential  equation  of  first 
order  for  J(s)  »  L[J(t) ]  with  variable  coefficients  of  not  more  than  third  order 
in  s.  The  general  solution  can  then  be  obtained  by  inversion,  I(t)  *  L  ^[J(s)]. 

It  should  be  noted  that  in  the  derivation  of  the  integral  equations  (25) 
and  (28),  the  current  I(t)  has  been  assumed  to  flow  in  the  +x  -  direction  (Fig.  3). 
Accordingly,  one  exspects  solutions  I(t)  <  0  for  "+"polarization  and  solutions 
I(t)  >  0  for  ,,-H  polarization. 
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ANALYTICAL  SOLUTION 


For  electrical  power  generation,  the  main  interest  is  in  shocked  solids 
with  sufficient  pressure  ionization  so  that  6  =  tq/t  <<  1.  For  this  reason, 
analytical  solutions  are  derived  for  the  case  of  small  space  charge  relaxation 
times  compared  with  the  polarization  relaxation  time  t.  The  solution  for 
large  space  charge  relaxation  times,  6  =  tq/t  >>  1,  can  also  be  given  in  closed 
form  — } 

For  fast  space  charge  relaxation,  <<  t,  the  initial-value  problem  (28) 

-  (29)  reduces  for  the  polarizations  P  in  (+)  and  opposite  (-)  to  the  direction 
of  shock  propagation  to  (6  <<  1) : 

+  (R  +  *)J(*)  -  ±6(1  -  e'*)  ,  (30) 

5  0 

J(-t-O)  =  0  .  (31) 

Eq.  (30)  consists  of  a  cumulative  voltage  term,  the  ohmic  voltage  loss,  and  a 
delayed  stress  polarization  voltage  source. 

By  means  of  the  transformation, 

J(*)  “  d£  .  ¥(t-0)  -  0  ,  <|»’(t«0)  -  0  ,  (32) 

the  integral  equation  problem  (30)  -  (31)  is  reduced  to  an  initial-value  problem 
for  a  first  order  differential  equation: 
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(33) 


(R  +  £  +  j  (R  +  a)<p  =  ±6(1  -  e_;t)e'(e/6);t 

<H-t=0)  =  0  .  (34) 

Eq.  (34)  implies  d\p(£-0)dt=0  by  Eq.  (33).  The  solutions  of  Eqs.  (33)  -  (34) 
are  for  the  ,,±H  polarization  cases 


-f-(a+R)  t  f(a+R)-l  ,  .  y 

*«)  =  ±6  (R  +  t)  /  (R  +  ty  (1  -  e  )e*1  '  '  <1*  ,  0  <  t  <  CL. 


(35) 


The  corresponding  cv rrent  solutions  J(-t)  for  the  "±"  polarization  cases  are  by 
Eq.  (32): 

t  y«f(a+R)-l  j  ,  ...  j 

t  /  (1  +  |)6  (1  -  e'*)e~(e/S)*<Lt' 

jw/6  - ±-^~- ;V-  *  - - 


6R 


(1  +  f>* 


f(a+R)+l 


-(e/6)* 


0  <  t  <  a  .  (36) 


These  solutions  can  be  extended  to  times  t  >  a  beyond  the  transit  time  t  *  a  if 
the  assumption  is  made  that  the  shock  wave  is  not  reflected  at  the  electrode  x  •  a. 

Eq.  (36)  indicates  that  the  integral  term  decreases  much  faster  with  in¬ 
creasing  resistance  R  than  the  ordinary  stress  polarization  term.  The  "+" 
polarization  current  is  -  by  magnitude  -  as  large  as  the  "-"polarization  current, 

|J(?)|  -  |j($)|  ,  0  <  t  <  a.  .  (37) 

Since  the  integral  term  is  dominant  by  magnitude,  the  solutions  of  Eq.  (36)  are 
J (t)  2*  0  for  the  ±  polarizations  and  times  t  <  a. 
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At  the  transit  time  t  -  <x,  the  integral  term  in  Eq.  (36)  vanishes,  and 
at  this  particular  instant  the  current  is  given  by 

J (t~a)  =  ±6(1  -  e~a)/(R  +  a)  .  (38) 

The  solution  in  Eq.  (36)  shows  that  the  (dimensional  analysis)  estimate 

—1  2 

I  -  APv/a  (Table  I)  is  correct  for  t  <  a/v  and  R  *  a  a/A  -  R,  since  I  **  (ca/R  )I 

a  o 

2 

*  (c, /e0) (R  /R)  APv/a  by  Eq.  (36).  For  power  applications  external  resistances 
x  z  a 

R  -  R  of  the  order  of  the  slab  resistance  are  of  main  interest, 
a 

Fig.  4  shows  computer  calculations  of  current  pulses  ("±"  polariza- 

tions)  versus  time  0  <  t  <  a  for  R  3  1  and  a  ■  1,2,  ...  10  (6 /e  *  10  )  based  on 

Eq.  (36).  The  pulse  height  | J (t) | increases  considerably  with  increasing 
a  *  a/vT,  i.e.,  if  the  slab  width  a  is  made  larger  and  larger  compared  with  the 
distance  vr  travelled  by  the  shock  wave  during  the  relaxation  time  t  (saturation 
for  a  >>  vx) . 

The  generated  currents  J(£)  are  large  for  R  $  1,  but  decrease  considerably 
(about  «  R  with  increasing  R  >  1.  The  steep  rise  of  the  current  J (£)  within  a 
time  t  %  10  \  i.e.,  t  %  t/10  is  remarkable.  Apparently,  the  complete  build-up 
of  the  electric  polarization  field  is  not  required  for  the  development  of  the 
maximum  current  flow  J (t)  if  t  <<  x  ,  since  the  electric  polarization  fields  are 
already  very  large  for  t  *  x/10. 
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Current  pulse  +  J (t)/6  (M±M  polarization)  versus  time  t  for 


CONCLUSION 


A  one-dimensional  theory  of  a  solid  state  power  generator  has  been  pre¬ 
sented,  in  which  the  electromotoric  force  is  produced  by  electric  polarization 
of  a  solid  in  the  stress  field  of  a  shock  wave.  The  numerical  results  indicate 
that  energies  of  the  order  AE  -  10^  Joule  can  possibly  be  generated  per  pulse. 
This  type  of  power  generator  does  not  require  an  external  magnetic  field  and 
appears  to  be  promising  as  a  pulsed  power  source  for  special  applications 
(for  which  the  efficiency  of  transformation  of  chemical  energy  of  explosives 
into  stress  waves  and  electrical  energy  is  of  no  concern) . 

Experiments  should  be  carried  through  to  verify  the  theoretical  expec¬ 
tations.  In  comparing  experimental  and  theoretical  current  pulse  forms,  it 
should  be  noted  that  the  temporal  current  distribution  depends  noticeably  on 
the  shape  of  the  stress  wave.  Current  pulse  forms  for  other  stress  wave  shapes 
than  the  simple  step  wave  can  be  calculated  in  complete  analogy.  It  might  also 
be  of  interest  to  refine  the  theory  by  including  multi-dimensional  effects  and 
deformation  of  the  solid. 
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ABSTRACT 

The  initial-boundary-value  problem  for  the  diffusion  of  an  in¬ 
itially  homogeneous  magnetic  field  into  a  slab  of  conductivity  a  <  00 
and  width  Ax  *  2a  is  solved,  under  consideration  of  the  electromagne¬ 
tic  wave  pulses  generated  at  the  surfaces  of  the  conductor  by  its 
interaction  with  the  external  magnetic  field,  which  propagate  into 
the  surrounding  vacuum.  The  analytical  solutions  show  that  (i)  the 
external  electromagnetic  transients  are  necessary  in  order  to  correc¬ 
tly  satisfy  the  boundary  conditions  for  the  tangential  electric  and 
magnetic  field  components,  and  (11)  the  spatial  and  temporal  develop¬ 
ment  of  the  electromagnetic  field  and  electric  current  in  the  conduc¬ 
tor  is  quantitatively  determined  by  a  new  dimensionless  parameter  group 

-1/2 

R  *  UQ<rac  [c  *  (uqco)  1*  This  "magnetic  Reynolds  number  of  the  va¬ 
cuum"  determines  the  coupling  between  the  transient  fields  in  the  con¬ 
ductor  o  >  0  and  the  ambient  space  (a  -  0). 


INTRODUCTION 


The  diffusion  of  electromagnetic  fields  B(r»t),  E(r,t)  in  a  conductor  of  finite 

conductivity  a  and  normal  surface  vector  n(s),  when  the  electromagnetic  field  ^(r,t) 

and  ^  (r,t)  outside  of  the  conductor  are  known,  is  in  general  described  by  Maxwell's 
o 

equations  without  displacement  current,  where  the  tangential  field  components  are 
assumed  to  satisfy  the  boundary  conditions'*^  n  *[B(r,t)  -  BQ(r,t)]  *  t)  and  S  * 

[£(r,t)  -l5o(r,t)]  *  t).  If  the  external  electromagnetic  field  is  time-independent 

and  electric  potential  fields  are  absent,  then  Bq  *  BQ(r)  ^  Eo  *  ^  (since  V  x  Eq  * 

^  ^  ^  12) 
-3B  /3t  *  0  and  E  *  r  0),  so  that  the  tangential  boundary  conditions  reduce  to  * 

o  o  o 

n  x  [B(r,t)  -  BQ(r)]  ■  t)  and  n  x  E(r,t)  *  ^ •  These  boundary  conditions  have  found  wide- 

1)  2) 

spread  use  in  mathematical  physics  ,  eleccromagnt  ic  theory  ,  and  the  theory  of  mag- 

3  4) 

netic  flux  compression  (at  the  outside  surface  of  the  liners)  *  .  However,  these 

boundary  conditions  are  questionable  approximations,  since  they  do  not  take  into  con- 

sideration  the  wave  fields  B(r,t),  E(r,t)  propagating  away  from  the  conductor  into  the 

-►  -i 

surrounding  medium,  which  have  their  sources  in  the  transient  current  fields  j  * 

V  x  b  of  the  conductor. 

For  a  concrete  illustration  of  the  problematics,  consider  the  diffusion  of  an 
external  (homogeneous)  magnetic  field,  Bq  *  {0,Bq,0}  for  |x|  >  a,  into  a  conducting 
slab  in  the  region  |x|  <  a  which  is  field  free  at  time  t  *  0  (Fig.  1).  Using  the  con¬ 
ventional  boundary  conditions,  the  transient  magnetic  field  B(x,t)  *  (0,B(x,t),0)  in 
the  conductor  is  determined  by  the  parabolic  initial-boundary-  value  problem:^ 

3B/3t  -  k32B/3x2  ,  |x|  <  a  ,  t  >  0  (1) 

B(x,t»0)  -  0  ,  |x|  <  a  ,  (2) 

B(x*±a,t)  *  B  ,  t  >  0  ,  (3) 

o 


where  <  *  1/uqo;  By  means  of  Fourier's  method,  the  general  solution  of  Eqs.  (1)  -  (3) 
Is  obtained  es  a  superposition  of  eigenfunctions^ : 


t  >  0,  (4) 


with  B(x,t)  Bq  in  |x|  <  a  for  t  ■+•  00 .  Since  V  x  B  *  u0aE»  electric  field 
E(x,t)  *  {0,0,E(x,t)>  in  the  conductor  is 


2B 


E(x,t)  = 


y 


,  iXn-l  -ic(2n-l)^Tr^t/4a^  u  2n-l 
(-1)  e  sin  m~2  - 


<  a,  t>  0. 


(5) 


In  accordance  with  the  boundary  conditions  (3),  the  space  outside  of  the  conduc¬ 
tor  remains  unperturbed  while  the  electromagnetic  field  diffuses  into  the  conductor, 

BQ(x,t)  *  Bq  ,  EQ(x,t)  -  0  ,  |x|  >  a  ,  t  >  0  (6) 

The  transient  currents  j  «  V  x  B/yQ  in  the  conductor  are  "eddy  currents and, 
therefore,  cannot  produce  transient  magnetic  fields  Bo(x,t)  -  BQ(x,t)  -  Bq  ^  0 
in  the  outside  region  |x|  >  a.  The  net  current  I(t)  through  any  cross  section 
2  *  constant  vanishes,  due  to  the  boundary  conditions  (3): 

1  B 

I (t) / Ay  *  J  *  [3B(x,t)/3x]dx  -  \T  /  °  dB  *  0  (7) 

-a  B 

o 


By  comparing  the  conductor  solutions  (4)  and  (5)  with  the  vacuum  solutions  (6), 

it  is  seen  that  B(x  *  ±a,t)  -  Bq  ~  0,  but  E(x  *  ±a,t)  -  Eq(x  *»  ±a,t)  *  E(x  ■  ±a,t)  #  0! 

1  2) 

Thus,  the  conventional  boundary  conditions  9  lead  to  a  violation  of  the  fundamental 
law  of  the  continuity  of  the  tangential  electric  fields  at  interfaces. 

The  correct  formulation  of  the  boundary  conditions  requires  consideration  of 
the  simultaneous  wave  fields  B+(x,t),  E+(x,t)  propagating  with  the  speed  of  light  c 
into  the  positive  and  negative  half-spaces  x  >  +a  and  x  <  -a  surrounding  the 
conductor  (Fig.  2),  which  are  excited  by  the  transient  current  fields  j(x,t)  « 
u~Hb(x, t)/9x  in  the  (space-charge  free)  conductor.  No  matter  how  small  the  external 
transients  B±(x,t)  and  E±(x,t)  are  (in  comparison  with  Bq  i*  0  and  Eq  *  0) ,  they  have 
to  be  taken  into  account  in  order  to  rigorously  satisfy  the  boundary  conditions 
n  *  [5]  ■  5  and  n  *  [E]  ■  5  for  the  continuity  of  the  tangential  electromagnetic 
fields  at  conductor  Interfaces. 
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The  quantitative  assessment  of  the  significance  of  the  external  transients  of 
the  diffusion  process  leads  to  the  discovery  of  a  new  dimensionless  parameter  combi¬ 
nation,  which  has  the  physical  meaning  of  a  "magnetic  Reynolds  number  of  the  vacuum": 

—  1/2  8 

R  *  UQcrac  ,  c  a  (uQeQ)  *  3  *  10  m/sec  .  (8) 

In  the  following,  the  formulation  of  the  initial-boundary-value  problem  for 

diffusion  processes  with  external  transients  and  its  analytical  solutions  for  the 

transient  electromagnetic  fields  inside  and  outside  the  conductor  are  presented.  The 

qualitative  and  quantitative  importance  of  the  new  boundary  conditions  and  the 

external  wave  fields  are  discussed  in  terms  of  R. 

The  presented  theory  has  important  implications  for  the  evaluation  of  the  flux 

3  4) 

losses  through  the  liners  of  magnetic  field  compressors  9  ,  the  electromagnetic 

acceleration  of  conducting  macroparticles^* the  electromagnetic  induction  in 

8  9) 

conductors  moving  relative  to  external  magnetic  fields  *  ,  and  for  the  interaction 

of  transient  plasma  shock  waves  with  external  magnetic  f ields10*11^ .  The  general 
significance  for  theoretical  physics  is  obvious. 
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INITIAL-BOUNDARY-VALUE  PROBLEM 


I 

I 

I 


I 


Subject  of  the  consideration  is  the  diffusion  of  the  magnetic  field  into  a  con¬ 
ducting  slab  |x|  <  a,  which  is  initially  embedded  in  a  homogeneous  magnetic  field 

B  *  (0,B  ,0,},  under  simultaneous  emission  of  electromagnetic  waves  from  the  conduc- 
o  o 

tor  surfaces  x  *  ±a  (Fig.  1). 

The  transient  electromagnetic  fields  =  {0,B+ (x,t) ,0}  and  H+  ■  {0,0, E+ (x,t)} 
in  the  infinite  vacuum  half -spaces  (o  =  0,  y  =  uq)  x  >  +a  and  x  <  -a  are  determined 


by  the  initial-boundary-value  problems  (±)  for  the  wave  equation, 

32B+/3t2  •  c232B+/3x2,  ±x  >  a,  t  >  0,  (9) 

B  (x,t  *  0)  *  B  ,  ±x  >  a  ,  (10) 

±  O 

B+(x  =■  ±a,t)  *  Bq  +  i|»+(t)  ,  t  >  0,  (11) 

since 

3E+/3 1  *  c23B+/3x  ,  3 E+  / 3 x  =  3B+/3t  (12) 

by  Maxwell's  equations  with  displacement  current.  The  solutions  of  Eqs.  (9)  -  (11) 
for  the  still  undetermined  boundary  values  i^+(t)  are: 


and 


B+(x,t)  -  B  +  *+(t  +  ~  ) 


-  B 


a  <  ±x  <  a  +  ct 
a  +  ct  <  ±x  <  ® 


E+(x,t)  «  +c<|/+(t  +  ) 

-  0 


a  <  ±x  <  a  +  ct 
a  +  ct  <  ±x  <  • 


(13) 


(14) 


These  solutions  are  typical  for  hyperbolic  equations,  i.e.,  the  boundary  values  <|<+(t) 
are  "transported"  into  the  half-spaces  ±x  >  a  with  the  speed  of  light  c,  so  chat  dis¬ 
continuous  wave  fronts  result  at  x  ■  ±(a  +  ct). 
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Let  the  external  magnetic  field  Bq  be  switched  on  at  t  3  0  so  that  no  electro- 

1  2) 

magnetic  fields  exist  in  the  conductor  for  t  <  0  9  .  The  conductor  has  a  finite  con- 

12)  ■+■  -4* 

ductivity  a  and  can,  therefore,  not  carry  surface  current  densities,  j*  *  ^ig^c EAx  =  0 
for  a  <  00  and  E  bounded.  Accordingly,  the  boundary  conditions  for  the  tangential  electric 
and  magnetic  field  components  at  the  conductor  vacuum  interfaces  are 


where 


B(x  -  ±a,t)  -  Bq  + 

,  t  >  0 

(15) 

E(x  =  ±a,t)  =  +ci|»+(t) 

,  t  >  0  , 

(16) 

E(x,t)  *  ic3B(x,t)/9x 

,  |x|  <  a  ,  t  >  0 

(17) 

by  Ohm’s  law  is  the  electric  field  in  the  conductor,  and  B(x,t)  is  the  magnetic  field 
in  |x|  <  a.  Furthermore 


k  -  1/u  a  >  0  .  (18) 

o 

(The  boundary  conditions  n«[cE]  =  0  and  n*  [B]  *  t)  are  satisfied  since  ^  and  1*  have  no 
normal  components.)  By  elimination  of  the  unknown  boundary  values  E(x  *  ±a,t)  and  \p+  (t) 
from  Eqs.  (15)  -  (17),  boundary  conditions  involving  only  the  magnetic  field  B(x,t)  in 
the  conductor  are  obtained: 


±a,t) 


t  >  0 


(19) 


These  are  the  new  boundary  conditions  for  the  diffusion  of  magnetic  fields  B(x,t)  into 

1  2) 

conductors.  They  differ  from  the  conventional  boundary  conditions  *  through  the  curl 

terms  3B(x  -  ±a,t)/3x  f  0,  which  consider  the  emission  of  magnetic  dilution  waves  from 

the  conductor  surfaces  x  *  ±a  into  the  vacuum  [ x|  >  a. 

Within  the  conducting  slab  of  finite  width  2a,  the  propagation  of  the  magnetic 

1  2  J2) 

field  can  be  treated  in  the  nonrelativistlc  or  diffusion  approximation  9  ^  .  Accor¬ 
dingly,  B(x,t)  in  the  initially  field  free  conductor  is  determined  by  the  parabolic 
initial-boundary-value  problem: 
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2  2 

3B/3t  *  Kd  B/3x  ,  |xj  <  a  ,  t  >  0  , 


B(x,t  -  0)  =  0 


9 B(x  =  ±a,t)/ax  ±  hB (x^+a, t)  =  ihB^,  t  >  0 


where 


h  ■  c/k  >  0 


The  transformation, 

B(x,t)  *  +  B(x,t)  ,  |  xj  <  a  ,  t  >  0  ,  (24) 

reduces  Eqs.  (21)  -  (22)  to  an  initial-boundary-value  problem  with  standard  "radiation" 
boundary  conditions: 

3B/3t  *  k3^B/3x^  ,  |x|  <  a  ,  t  >  0  ,  (25) 


B(x,t  *  0)  *  -B  ,  |x|  <  a 


3S(x  =  ±a,t)/ax  ±  hB(x=s±a, t)  =»  0 


In  accordance  with  Fourier's  theorem,  the  solution  of  Eqs,  (25)  -  (27)  Is  obtained 
as  a  superposition  of  eigen-functions  B^(x,t)  for  the  region  j x|  <  a  which  satisfy  the 
boundary  conditions  (27): 

*  (h^+k^)k  *sink  a  -<k2t 

B(x,t)  *  -2B  l  - - — ? - —  e  n  cosk  x  ,  |x|  <  a  ,  t  >  0  ,  (28) 

n-1  [Q»-h£)*Hi]  n 

where 

k^a  tg(kfta)  *  ha  ,  n  -  1,2,3*.*.,  (29) 

gives  the  eigenvalues  ^associated  with  the  boundary  conditions  (27), 

1,  Conductor  Solutions: 


For  a  compact  representation  of  the  field  solutions,  dimensionless  independent 


(£,t)  and  dependent  variables  are  Introduced, 
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5  =  x/a 


(30) 


=  Kt/a 


a  =  k  a 
n  n 


8(5, x)  =  B(x,t)/BQ,  E(5,x)  =  E(x,t)/(i<Bo/a) ,  =  j (x,t)/(BQ/yoa). 

(31) 

According  to  Eqs. (24)  and  (28),  the  solutions  for  the  dimensionless  fields  8(5  ,t), 
E(5,x)  *  38(5»x)/35,  and  /(5,x)  in  the  conductor  are  given  by: 

,p2.  2.  -1  .  2 
00  (.k  +a  )a  sina  -a  x 

8(5, t)  -  1  -  2  J  - f-f - —  e  n  cosaj  ,  |?|  <  1  ,  t>  0  ,  (32) 

n=l  [  (R  -hx^)+R]  n 

,„2.  2.  .  2 
“  (k  +a  )  sina  -a  T 

E( 5,t)  =  2  £  - 5^-s - —  e  n  sina  5  ,  |5|  <  1  ,  t  >  0  ,  (33) 

n-1  [(RZV)+R]  n 

j(5» t)  -  E(5,x)  ,  |cf  1  1  ,  t  >  0  ,  (34) 

where 


a  tga  =  R  ,  n  =  1,2,3,...,  R  *  ac/<  =  ti  oac 
n  n  o 


(35) 


-2 

by  Eqs.  (8)  and  (29).  For  sufficiently  large  times  t  >>  ,  the  homogeneous  magnetic 

field  has  diffused  completely  into  the  conductor. 


8(£,t)  -►  1  ,  E(£,t)  +  0  ,  /(£, t)  -►  0  ,  t  +  » 


(36) 


In  the  hypothetical  limit  of  infinite  magnetic  vacuum  Reynolds  number  R,  Eq.  (32)  re¬ 
duces  to  the  known  solution  (4)  for  the  conventional  boundary  conditions'^ , 


00 

limB(5  ,t)  -  1  *  7  I  -^TT 
R-~  irn«l',2n  1} 


«  ,  ,,n-l  _(2n-l)2iT2T/4  2n-l  _  ^  ^  A 

e  cos — —  ir£,  |5|  <  1,  t  >  0, 


since 


lim  a  *  2°.^  tt,  n  ■  1,2,3,... 
n  2 


(37) 

(38) 


Comparison  of  Eq.  (32)  with  Eq.  (37)  indicates  that  the  8(5  ,x)  solutions 


with  the  new  and  conventional  boundary  conditions  differ  not  much  if  R  >>  a^  ■  iv/2. 
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2.  Vacuum  Solutions: 


In  view  of  the  boundary  conditions  (15),  Eq.  (32)  yields  for  the  boundary 
values  ¥+(t)  =  8(£=±1,t)  -  1.  Accordingly,  Eqs.  (13)  and  (14)  give  for  the  electro¬ 
magnetic  fields  in  the  positive  (£  >  +1)  and  negative  (£  <  -1)  half-spaces  the  so¬ 
lutions  ; 


and 


8±(C,t)  =  1  +  Y±(t  +  ,  1  <  ±£  <  1  +  Rt  , 

=  1  ,  1  4*  Rt  <  ±£  <  °°  , 


E±(e, t)  =  +RV(x  +  ,  1  <  ±5  <  1  +  Rt  , 


=  0 


1  +  Rt  <  ±£  < 


where 


(t  +  -  -2  J 


(R2+a^)sinancosan  -aJ[T+<€+l)/R] 


n=l  an[(R2+a^)+Rl 


(39) 


(40) 


(41) 


Eqs.  (39)  and  (40)  represent  electromagnetic  wave  pulses  which  propagate  with  the 
dimensionless  speed  R(c)  from  the  conductor  surfaces  £  -  ±1  into  the  vacuum  spaces 
±£  >  1  with  discontinuous  wave  fronts  at  £  =  ± (1  +  Rt).  They  are  kicked  on  at  t  3  0 
and  their  emission  lasts  to  the  end  (t  -►  «)  of  the  diffusion  process  in  the  conductor. 
The  vacuum  fields  8+(£,t)  are  in  opposite  direction  of  Bq ,  i.e.,  they  represent  di¬ 
lution  waves  (Fig.  2). 

In  the  case  of  large  coupling  numbers,  R  >>  1,  Eqs.  (39)  and  (40)  yield  with 


Eq.  (41) 


-1, 


8+(£,t)  *  1  +  0[R 


E+  (£  ,t)  *±2  f  e 
n-1 


,2n-l  .2 
-C  2 — B  )  x 


1  <  ±C  <  1  +  Rt  ,  t  >  0,  R  »  1,  (42) 
1  <  ±5  <  1  +  Rt  ,  t  >  0,  R  >>  1,  (43) 


since  cosa  * 
n 


(-l)n_1(2n-l)W2R  for  R  »nuby  Eq.  (35). 
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The  magnetic  field  Bq  outside  of  the  conductor  remains  nearly  unperturbed 
during  the  diffusion  process,  -  R  \  whereas  the  external  electric  transients 
^  0  are  of  order  behind  the  wave  fronts,  £  =  ±(1  +  Rt),  for  R  >>  1.  However, 
since  V*B+  38  c  ^3E^/3t,  not  only  but  also  -  R  ^  cannot  be  neglected  for 

R  >>  1. 

Thus,  we  have  shown  how  self-consistent  solutions  can  be  obtained  for  the 

electric  and  magnetic  fields  in  the  conductor  and  the  surrounding  vacuum,  which 

satisfy  the  boundary  conditions  for  the  continuity  of  the  tangential  electric  and 

magnetic  fields  at  the  conductor-vacuum  interfaces.  The  conventional  boundary 

1  2) 

conditions  for  electromagnetic  diffusion  processes  *  ,  ignore  the  external  elec¬ 

tromagnetic  transients,  violate  the  boundary  condition  for  the  tangential  electric 

field,  and  permit  no  Poynting  vector  S  =  ExH  outside  of  the  conductor.  As  a  result, 

1  2) 

the  conventional  boundary  conditions  9  make  it  impossible  for  electromagnetic 

fields  to  diffuse  through  condr.rc.  r  and  to  escape  into  the  ambient  space. 

For  both  the  conductor  an  solutions,  the  limit  R  0,  which  implies 

a  -*  0  since  a  t  0,  is  not  realizable  since  the  conductivity  of  conductors  is  by 

definition  large.  For  insulators  or  extremely  poor  conductors  (R  *  a  ■+-  0) ,  the 

12) 

nonrelativistic  or  parabolic  diffusion  equation  is  inapplicable  .  Therefore, 
the  investigation  of  the  limit  R  ->■  0  would  require  solution  of  Maxwell’s  equations 
with  displacement  current  in  the  slab  |x|  <  a  given  elsewhere^^ . 

The  generalization  of  the  theory  required  for  conductors  and  external  media 
(vacuum,  gases,  fluids)  with  different  permittivities  e  and  u  are  trivial  but 
complicate  the  notation. 
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DISCUSSION 


where  £  is  the  extension  of  the  conductor.  For  conductors,  the  diffusion  time  is 

-7  7  -4 

relatively  large  if  £  is  not  microscopically  small,  e.g.,  *  4tt*10  6*10  *10 

-2  -2 
~  10  sec  for  a  copper  slab  of  width  a  *  10  m. 

2+  2 

Comparison  of  the  neglected  term  3  B/3t  with  the  second  and  third  terms  of 
Eq.  (44)  reveals  the  relation  of  the  parabolic  diffusion  approximation  to  the  new 
coupling  number  R  *  p  aac: 


6b 


-1  3Bi 
R  St1 


XR  Vo  p-2 
—  *  - r  *  k 

TD  (u  aa) 
o 


(48) 


This  result  again  confirms  the  validity  of  the  parabolic  diffusion  equation  for 

conductors,  for  which  R  =  U^aac  >>  1.  E.g.,  R  =  4ttx10~76x107x10~23x108  -  108 

-2 

for  a  copper  slab  a  3  10  m.  More  important,  Eq.  (4g)  demonstrates  that  the 

2+  2 

neglected  relativistic  term  3  B/3t  in  the  conductor  is  small  of  order  R  <<<<  1, 
whereas  the  calculated  electromagnetic  fields  in  the  conductor  are  of  order 
8  -  E  -  R^  [Eqs.  (32)  -  (33)],  and  the  external  electromagnetic  transients  are 
of  order  B+  -  R  ^  and  E+  -  R^  [Eqs.  (39)  -  (40)],  since  in  Eq.  (41)  for  large  R 

|  cosaj  *  [1  +(R/a^)^]  *  cWR  ~  ^  ,  R  »  mr  ,  n  ■  1,2,3,...  . 

(49) 

In  conclusion,  it  is  noted  that,  in  conductors,  magnetic  field  diffusion  is 
a  nonrelativistic  process  (as  is  electric  conduction,  j  *  oE) .  The  electric  tran¬ 
sients  E+  in  vacuum  must  be  of  the  same  order  as  the  electric  field  E  in  the  con¬ 
ductors,  E+  -  E  -  R°,  since  otherwise  the  tangential  electric  field  could  not  be 

continuous  across  the  conduct or-vacuum  interface.  On  the  other  hand,  the  external 

~  -1  —1  —1 
magnetic  transients  B+  are  small  of  order  R  *  (pQoa)  c  since  the  magnetic 

field  energy  flows  with  the  speed  of  light  towards  the  conducting  cavity.  The 

deeper  physical  reason  for  these  electromagnetic  transients  is  to  be  seen  in  the 

conservation  laws  for  electromagnetic  energy  and  momentum,  which  follow  from  Max- 

12) 

well’s  equations 
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INITIAL-BOUNDARY -VALUE  PROBLEM  FOR  ELECTROMAGNETIC  INDUCTION  IN  ACCELERATED  CONDUCTORS 

MOVING  ACROSS  MAGNETIC  FIELDS 


H.  E.  Wilhelm 
Michelson  Laboratory 
Naval  Weapons  Center 
China  Lake,  California  93555 


Abstract 

Boundary  conditions  are  derived  for  the  interfaces  of  a  conductor  moving  across 
an  external  magnetic  field  in  an  ambient  medium  (vacuum  or  non-conductor),  which 
consider  the  emission  of  electromagnetic  waves  from  the  conductor  surface  as  a 
result  of  electromagnetic  induction.  These  boundary  conditions  are  applied  to 
the  initial-boundary-value  problem  for  the  electromagnetic  induction  in  a  con¬ 
ducting  slab,  which  is  accelerated  across  a  homogeneous  magnetic  field  to  a  non- 
relativistic  velocity.  Fourier  series  solutions  are  presented  for  the  transient 
electromagnetic  fields  in  the  moving  conductor  and  the  discontinuous  electromag¬ 
netic  waves  in  the  ambient  space.  It  is  shown  that  the  transient  electromagnetic 
fields  inside  and  outside  the  conductor  are  due  to  two  mechanisms,  i.e.,  "velocity 
induction"  (ordinary  induction)  and  "acceleration  induction"  (dv(t)  / dti^O ) . 
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INTRODUCTION 


Although  the  theoretical  foundations  for  electromagnetic  induction  in  con¬ 
ductors  moving  across  magnetic  fields  were  formulated  in  1908  by  H.  Minkowski^, 

only  simple  problems  such  as  stationary  unipolar  induction  in  rotating  discs 

2) 

have  been  discussed*  The  electrodynamics  of  moving  media  has  been  a  subject 

3) 

of  basic  research  and  also  controversy  to  date.  In  the  treatment  of  electro- 

4  5) 

magnetic  induction,  e.g.,  in  liners  of  magnetic  field  compressors  *  and  tran¬ 
sient  plasma  shock  waves  interacting  with  external  magnetic  fields  B  it 

has  become  customary  to  use  the  boundary  condition  for  the  tangential  magnetic 
field  in  the  form  nx[  B  -  Bq]  =  at  the  conductor-gas  interface,  where  B  is 
the  transient  magnetic  field  in  the  conductor  and  Bq  is  the  unperturbed  (!) 
external  magnetic  field.  This  boundary  condition  leads  to  electromagnetic 
solutions  B,E  in  the  conductor  which  do  not  satisfy  the  corresponding  bound¬ 
ary  condition  uxf  E  -  ?q]  •  (J  for  the  tangential  electric  fields  at  the  mov¬ 
ing  interface.  These  "conventional"  boundary  conditions  produce  approximate 
to  incorrect  results,  depending  on  the  physical  situation. 

For  an  illustration  of  the  problematics  of  the  conventional  boundary 
4  5) 

conditions,  ’  which  are  also  being  used  in  the  analysis  of  magnetic  field 

8) 

diffusion  into  conductors  at  rest  ,  consider  a  conducting  slab  Ax  -  2a  with 
its  surfaces  initially  at  x  -  ±a  in  a  transverse  homogeneous  magnetic  field 
?o  *  {0,Bq,0}  for  |x|  <  »  (Fig.  1).  At  time  t  -  0,  this  conductor  is  set 
into  motion  with  a  velocity  v(t)  *  (d£(t) /dt, 0,0}  so  that  its  front  and  rear 
surfaces  are  at  x  *  x(t)  t  a  at  time  t  >  0  where  £(t-0)  ■  0.  No  matter  whether 
the  induction  of  the  transient  magnetic  field  ?(x,t)  in  the  moving  conductor 
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■Hr-; 


9) 

is  described  by  the  relativistic  wave  equation  or  the  non-relativist ic  dif- 
fusion  equation  ,  the  initial  condition  B(x,t=0)  =  B^,  -a  <  x  <  +a,  and  the 
conventional  boundary  conditions  B(x=x(t)±a, t)  *  Bq,  t  >0,  permit  only  one 
and  the  same  solution,  B(x,t)  *  Bq,  which  impliej  E(x,t)  *  -v(t)xBQ,  and 
T(x,t)  *  (J  by  Ohm’s  law  for  the  moving  conductor.  These  simple  solutions  are 
due  to  the  conventional  boundary  conditions  without  external  perturbations  and 
are  obviously  not  correct,  since  the  boundary  condition  for  the  tangential  elec¬ 
tric  field  is  not  satisfied,  nx[  ?  -  E  ]  3  v(t)B  e  ^  3*  where  ?  *  S.  Since  M. 

o  o  y  o 

Faraday  it  is  an  experimental  fact  that  transient  electromagnetic  fields  and 
currents  are  induced  in  the  conducting  slab  as  soon  as  it  is  moved  relative  to 
the  external  magnetic  field  Bq. 

In  the  following,  an  analytical  solution  is  presented  for  the  initial- 

boundary-value  problem  of  the  electromagnetic  induction  in  a  conducting  slab 

Ax  =«  2a,  which  is  at  rest  for  t  <  0  in  a  transverse,  homogeneous  magnetic 

field  B  ,  and  which  is  set  in  motion  at  t  -  0  with  a  nonrelativist ic  velocity 
o’ 

v(t)  of  arbitrary  (finite)  acceleration  dv(t)/dt  (Figs.  1,2).  The  electro¬ 
magnetic  induction  In  the  conductor  is  shown  to  produce  transient  electro¬ 
magnetic  fields  within  the  conductor  and  electromagnetic  waves  at  the  moving 
conductor  surfaces  x  -  £(t)  ±  a,  which  propagate  with  the  speed  of  light  in 
the  surrounding  space  to  infinity.  Only  if  the  electromagnetic  waves  outside 
of  the  moving  conductor  are  taken  into  consideration,  consistent  solutions  of 
Maxwell's  equations  exist  which  satisfy  the  boundary  conditions  for  the  tan¬ 
gential  and  normal  electromagnetic  fields  at  the  conductor-non-conductor 
interfaces . 

The  analytical  solutions  for  the  moving  conductor  and  surrounding  spaces 

permit  significant  conclusions.  The  transient  electric  field  induced  in  the 
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conductor  is  the  sum  of  a  field  which  is  proportional  to  the  velocity  v(t) 
and  a  field  which  is  an  integral  functional  of  the  acceleration  dv(t)/dt 
of  the  conductor,  A  fundamental  dimensionless  group  (a  *  conductivity, 
permeability  of  conductor) 


R  =»  u-joac,  c 


^2e2^ 


-1/2 


is  found  which  represents  a  "magnetic  Reynolds  number"  of  (i)  "free  space" 

(if  the  conductor  moves  in  vacuum  or  gas  with  permeabilities  e0>^0,)»  or  (ii) 
"non-conducting  space"  (if  the  conductor  moves  in  a  nonconducting  medium  with 
permeabilities  R  determines  the  coupling  of  the  transient  electromagnetic 

fields  inside  and  outside  the  conductor.  The  external  magnetic  transients 
would  be  negligible  for  R  -►  00 ,  however,  the  external  electric  transients  do 
not  vanish  for  R  «.  Thus,  the  conventional  boundary  conditions  which  assume 
that  the  external  magnetic  field  Bq  remains  unperturbed  outside  of  the  conductor 
while  transient  electromagnetic  processes  occur  in  the  latter  are  not  correct 
within  electromagnetic  theory. 
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BOUNDARY  CONDITIONS 


For  the  analysis  of  electromagnetic  induction  in  moving  conductors  (conducti¬ 
vity  a  <  00 ,  permittivities  which  move  with  a  velocity  field  v(r,t)  relative 

to  the  "laboratory  system"  in  a  non-conducting  medium  (vacuum,  gas,  or  fluid;  o  *  0, 
z29^2^  9  t*lfe  boundary  conditions  for  the  electromagnetic  fields  ^  (conductor) 
and  ^2  (non“conductor)  are  required  in  the  L-frame.  Integration  of  Maxwell’s 
equations  with  displacement  current over  the  interface  1-2  with  velocity  v(r,t) 
and  normal  vector  n(r,t)  (direction  1  2)  yields  for  the  tangential  and  normal  elec¬ 

tromagnetic  field  components  the  boundary  conditions  in  the  L-system: 


where 


nx[^2  -  =  +(n*v)[^2  - 

nx[^2/v2  -  -  (n.v)[e2^2  -  e^l 

”'^e2^2  “  el^l^  *  p 
n*  [^2  -  1^1  *  0 

-fie  *-*  if 

j  :  -  p  v  ,  p  lim  pAs 

A  9^0 


(1) 

(2) 

(3) 

(4) 

(5) 


are  the  surface  current  and  charge  densities  of  the  Interface  As  0.  A  conductive 

surface  current  density  does  not  exist  at  a  conductor  of  finite  conductivity,  "f  * 

lim  allAs  ■  15  for  a  <  *>  and  1i  bounded* 

As+0 

In  many  cases,  the  permittivities  of  good  conductors  and  their  ambient  atmos¬ 
pheres  equal  the  free  space  values,  2  *  e  V1  2  *  With  this  simplification 
in  notation,  the  boundary  conditions  (1)  -  (4)  are  applied  to  the  front  (+)  and  rear 
(-)  surfaces  x  *  x(t)  ±  a  of  a  conducting  slab  (o  <  ®)  with  the  fields  (Fig*  1) 
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B  ■  {0,B(x,tXO)  ,  E  -  {0,0,E(x,t)}  ,  x(t)  -  a  <  x  <  x(t)+a  ; 

B(x,t=*0)  =  Bq  ,  E(x,t=0)  -  0  ,  -a  <  x  <  +a,  ,  (6) 

which  moves  with  the  velocity  v(t)  •  { dx(t)/dt,0,0}  in  an  ambient  medium  ( a  m  0)  with 
the  fields  (Fig.  1) 

B+  *  {0,B±(x,t),0}  ,  E±  *  {0,0,E±  (x,t)},  x  %  x(t)  ±  a  ; 

B+(x,t*0)  »  Bq  ,  l+(x,t-0)  =  3  ,  x  £  ±a  .  (7) 

BQis  an  external  homogeneous  magnetic  field  which  fills  uniformly  the  conductor  (1) 
and  the  medium  (2).  The  boundary  conditions  (3)  and  (4)  are  satisfied  identically 
since  B,E  and  B+,E+  are  1  n  so  that  "j  *  t)  and  p  *  0  by  Eq.  (5).  The  tangential 
boundary  conditions  (1)  and  (2)  yield  for  the  fields  (6)  -  (7) : 

E+(x*x(t)±a,t)  -  E(x*x(t)±a,t)  «  +v(t) [B± (x=x(t)±a,t)  -  B(x*x(t)±a,t)] , 

"  (8) 

B+(x»x(t)±a,t)  -  B(x-x(t)±a,t)  *  +  v(t)c~2[E+ (x»x(t)±a,t)  -  E(x-x(t)±a,t) ] , 

(9) 

where 

c  -  (ue)~1/2  (10) 

is  the  speed  of  light.  For  non-relativistic  conductor  motions*  Eqs.  (8)  -  (9)  reduce 
to: 

E+(x-x(t)±a,t)  -  E(x-x(t)±a,t)  -  0  ,  v(t)2  «  c2  ,  (11) 

B±(x-x(t)±a,t)  -  B(x-i(t)±a,t)  -  0  ,  v(t)2  «  c2  .  (12) 

According  to  the  non-relativistic  Ohm’s  law  for  conductors  which  are  moving  with 
a  velocity  v  relative  to  the  L-system*  1  *  a(1i  +  vx$)  *  the  electric  field  l£(x,t)  in 
the  conductor  is  expressed  in  terms  of  3(x,t), 
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—1  9 

E  *  -vB  +  (pa)  x3B/3x  ,  x(t)  -  a  <  x  <  x(t)  +  a  ,  v(tj 


(13) 


The  electromagnetic  field  (x,t) ,  E+ (x,  t)  in  the  adjacent  semi-infinite  half  spaces 
is  determined  by  the  hyperbolic  initial -boundary-value  problems: 

32B+/3t2  -  c232B+/3x2  ,  x  ^  x(t)  +  a  ,  t  >  0  ,  (14) 


B+(x,t=0)  -  Bq 


»  x  %  ±a 


since 


B+ (x*x(t)±a,t)  -  B(x»x(t)±a,  t)  , 


t  >  0  , 


3E+/3x  =*  3B+/3t 


,  3B+/3x  -  c  3E+/3t 


by  Maxwell’s  equations  for  homogeneous  non-conductors.  Eq.  (16)  couples  the  solu¬ 
tions  B+(x,t)  in  the  semi-infinite  spaces  x  *  x(t)  ±  a  to  the  solution  B(x,t)  in  the 
conductor,  x(t)  -  a  <  x  <  x(t)  +  a.  By  Eqs.  (14)  -  (17),  the  ambient  electromagnetic 
field  transients  are  of  the  form: 


x  +  a 


B,  (x ,  t )  -  B  +  ¥  (t  X - )  ,  x(t)  +  a  ^  x  i*  ±  (a  +  ct)  , 

+  o  ±  ^  c 


,  a  +  ct  <  ±x  <  » 


E+(x,t)  -  +c*+(t  +  .  x(t)  ±  a  *  x  9  ±(a  +  ct)  , 


where 


,  a  +  ct  <  ±x  <  ® 


*+< t  +  ^i)  -  B(x-x(t)±a,t)  -  Bfl 


by  Eq.  (16)  determines  the  form  of  the  wave  functions  f  (h  )  of  the  "self  similar" 
arguments  n  *  t  +  x/c  +  a/c  from  the  conductor  solution  B(x,t). 

The  solutions  (18)  -  (19)  are  typical  for  hyperbolic  initial-boundary-value 
problems,  i.e.,  the  boundary  values  ¥  (t  +  x(t)/c)  are  transported  with  the  speed  of 


light  c  into  the  half  spaces  x  ^  x(t)  ±  a,  so  that  discontinuous  wave  fronts  result 
at  x  *  ±(a  +  ct) . 

By  means  of  Eqs.  (18)  -  (20),  the  boundary-values  B+ (x=x(t)±a,t) ,  E+ (x=x(t)±a,  t) , 
and  (t  +  x(t)/c)  are  eliminated  from  the  boundary  conditions  (11)  and  (12): 

1  tc  +  v(t)]B(x*x(t)±a,t)  =  ±cBq,  v(t)2  «  c2  .  (21) 

This  are  the  fundamental  new  boundary  conditions  for  moving  conductors  which  (i)  in¬ 
volve  only  boundary  values  of  the  magnetic  field  B(x,t)  of  the  conductor,  and  (ii) 
consider  the  emission  of  electromagnetic  waves, 

B±(x,t)  *  ¥±<t  +  5~2-)  ,  £±(x,t)  *  +c«F(t  +  Z-~),  (22) 

from  the  conductor  surfaces  x  *  x(t)  ±  a  into  the  ambient  spaces  x  ^  x(t)  ±  a. 

Since  the  magnetic  field  in  the  conductor  is  the  sum  of  the  external  Bq  and  a  tran¬ 
sient  B(x,t), 

B(x,t)  -  Bq  +  B(x,t)  ,  x(t)  -  a  <  x  <  x(t)  +  a  ,  (23) 

Eq.  (21)  gives  for  the  transient  conductor  field  the  boundary  conditions: 

"~~k£ti£x±aj~?  i  +  v(t)]B(x»x(t)±a,t)  -  v(t)BQ,  v(t)2  «  C2  .  (24) 

4  -1 

In  view  of  the  difficulties  to  accelerate  macroscopic  bodies  to  speeds  )v|  >  10  m  sec 
the  stronger  non-relativistic  condition  |v|  «  c  is  in  general  fulfilled,  which  redu¬ 
ces  Eqs.  (21)  and  (24)  to  the  simpler  boundary  conditions: 

(po)-l  3B<Xr\^a,t?  ±  c[B(x*x(t)±a,t)  -  Bo]  -  v(t)Bo  ,  |v(t)|  «  c  , 


and 


j 

(uo)  1  ±  cB(x=x(t)±a,t)  -  v(t)B^  ,  |v(t)|  «  c  .  (26) 

If  the  Vxli  and  v(t)B^  terms  are  omitted,  Eqs.  (25)  and  (26)  reduce  to  the  con¬ 
ventional  boundary  conditions^  B(x=x(t)±a,  t)  *  Bq  and  B(x-x(t)±a, t)  *  0,  respec¬ 
tively.  Comparison  shows  that  the  conductor  currents  VxB/p  and  the  induced  currents 
crv(t)BQ  at  the  conductor  surface  are  the  sources  of  the  emitted  electromagnetic 
transients. 


IK IT I AL -BOUNDARY -VALUE  PROBLEM 


Consider  a  rigid  conducting  slab  of  width  &x  *  2a  with  surfaces  in  the  planes 
x  =  x(t)  ±  a  at  time  t  >  0  (Fig.  1).  This  conductor  is  exposed  to  an  external  mag¬ 
netic  field  *  {0»B^,0}  which  is  homogeneous  throughout  the  space  -<»  <  x  <  +°°,  and 
is  accelerated  to  a  (non- relativistic)  velocity  v(t)  -  {v(t),0,0}  from  an  initial  po¬ 
sition  at  rest  xCt^O)  5  0  so  that 


with 


v (t )  »  dx(t)/dt  ,  x(t)  *  /  v(t 1  )dt 1  ,  dv(t)/dt  £  0 

0 


x(t»0)  =  0  ,  v(t=0)  *  0  ,  dv(t=0)/dt  \  0 


,  (27) 

.  (28) 


The  electromagnetic  induction  of  the  transient  electromagnetic  fields  $  *  {0,B(x, t) ,0} , 
^  *  {0,0,E(x,t)>  in  the  conductor  of  finite  width,  as  a  result  of  its  accelerated  mo¬ 
tion  v(t)  across  the  external  magnetic  field  is  determined  by  the  parabolic  ini¬ 
tial-boundary-value  problem  for  B(x,t): 


If  +  v^ff  =  ~  a  <  x  <  *U)  +  a  ,  t  >  0 

3t  3x  ax2 


B(x,t-0)  -  B 


-a  <  x  <  +a 


,  (29) 

,  (30) 


t  hlB(x^(t)±a,t)  -  BJ  -  <_1v(t)B  ,  t  >  0  .  (31) 

oX  O  O 


where 


1/uo  ,  h  *  c/tc 


•  (32) 


Eqr.  (29)  follows  from  Maxwell's  equations  and  0hmrs  law  j  *  a(E  +  v*B)  for  moving  con- 

9) 

ductors  in  the  diffusion  approximation,  *  e/a  0.  The  boundary  conditions  (31) 
couple  the  electromagnetic  induction  process  in  the  conductor  to  the  external  transi¬ 
ents  in  the  ambient  medium  [Eq*  (25)]*  The  transformation. 


B(x,t)  -  Bq  +  B(X,t)  , 


-a  <  x  <  +a 
81 


t  >  0 


(33) 


t  >  0 


(34) 
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since  v(t*0)  *  0.  With  3F(x,t)/3t  odd  in  x  and  F(x,t*0)  -  0,  the  initial-boundary-va¬ 
lue  problem  (41)  -  (43)  is  solved  by  means  of  the  Fourier  expansions: 


vhere 


G(x,t)  3  I  G  (t)sin  k  x  >  -a  <  x  <  +a 
u  _  n  n 

n*l 


n*l 


k  a  ctg  k  a  »  -ha  ,  n  *  1,2,3,  ... 
n  n 


yields  by  means  of  the  inverse  Fourier  theorem 


dG  (t)/dt  +  Kk^G  (t)  -  -S  (t) 
n  n  n  n 


G  (t»0)  -  0 
n 

where 


Sn(t)  *  (a/K)KnBQdv(t)/dt 


K  -  (h2  +  k2)sin  k  a/(k  a)2[(h2  +  k2)  +  (h/a)] 
n  n  n  n  n 


t  -tck2 (t-t') 

G  (t)  -  -  /  e  n  S  (t')dt' 

n  0  n 


t  >  0 

» 

(47) 

t  >  0 

9 

(48) 

f 

(49) 

t  X  associated  with 

Che 

n 

)  into  Eqs.  (41) 

and 

l  (42) 

t  >  0 

* 

(50) 

(51) 

t  >  0 

* 

(52) 

a 

(53) 

is 

t  >  0 

• 

(54) 

Combining  of  Eqs.  (47),  (52),  and  (54)  yields  as  solution  of  the  initial-boundary-va¬ 
lue  problem  <41)  -  (43) 


2 

®  t  .  .  ,v  -*rk  (t-tf) 

G(x,t)  -  -  V  l  K  if  e  n  dt?]sin  k  X,  -a  <  X  <  +a,  t  >  0  . 

K  °n-l  n  0  dt  n 

(55) 

By  Eqs.  (33),  (34),  (38),  (44),  and  (55),  Che  magnetic  field  in  Che  moving  conductor 
is 
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BOc.t)  -  B  +  iMt)B l /2<Jc+ac)](x-x(t))  -  — B  f  K  { 

°  °  k  o  L  nvA  dt f  dt 

n*l  0 

x(t)  -  a  <  x  <  x(t)  +  a  ,  t  >  0 


]sinkn(x-x(t)) f 


(56) 


This  is  a  fundamental  result,  which  shows  that  the  transient  magnetic  field  B(x,t)  is 
the  sum  of  a  field  F(x,t)  induced  by  the  motion  v(t)  and  a  field  G(x,t)  induced  by  the 
acceleration  dv(t)/dt  of  the  conductor  in  the  external  magnetic  field  Bq.  Similar  de¬ 
compositions  exist  for  the  electric  field  E(x,t)  and  current  density  j(x,t)  in  the  ac¬ 
celerated  conductor. 

The  still  unknown  forms  (±)  of  the  two  external  wave  functions  ^+(t  +  x/c  +  a/c) 
in  the  spaces  x  ^  x(t)  ±  a  are  determined  from  the  solution  (56)  by  means  of  the  boun¬ 
dary  condition  (20) ,  which  gives 

*+(t  +  ^-)  *  B(x»±a,t)  ,  t  >  0  .  (57) 

The  transformations  t^  *  t+(t)  and  its  inverses  t  *  f+{t+)  for  the  two  (±)  waves  defined 
by 

,  t  -  f+{t±}  ,  t  >  0  ,  (58) 

where  t+  «  t_  *  0  for  t  *  0  (x(t*0)  *  0)  but  t+  i*  t  for  t  >  0,  show  that  the  wave  func¬ 
tions  are  of  the  fora  (t+)  »B()t*±  a,t  *  f+{t+>).  Accordingly, 

f+  (t  +  X-^-a)  -  B(X«±a,  t-f+(  t  +  ) ,  x(t)  ±  a  5  x  *  ±  (a  +  ct)  ,  t  >  0. 

(59) 

Since  B(x,t)  *  B(x,t)  -  Bq,  substitution  of  Eq. (56)  into  Eq.  (59)  gives  the  wave  func¬ 
tions  f  as  functionals  of  f+{t  +  x/c  +  a/c)  : 
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I 

I 


y±(t  +  “  4taBo/2(K  +  ac)]v(f±{t  +  S-J-*}) 

,2.  f  .  xf  a1  c  ^  .  2  , 

n  ±  x  c  J  dv(t')  n  ^  , 

i  dt’ 


x  —  B  F  K  sin  k  a  e 
T  k  o  n  it 
n*l 


e  “  dt'  , 


x(t)  ±  a  ^  x  &  ±  (a  +  ct)  >  t  >  0 


•  (60) 


The  solutions  (60)  determine  the  propagation  of  the  emitted  electromagnetic  waves 
outside  of  the  moving  conductor,  x  ^  x(t)  ±  a.  Again,  these  waves  consist  each  of  a 
"velocity"  wave  [v(t)J  and  an  "acceleration"  wave  [dv(t)/dt].  They  satisfy  all  boun¬ 
dary  and  initial  conditions, 

y±/x  -  x(t)  ±  a  =  *  t  >  0  ,  (61) 

Vt  »  o  *  0  ,  *  *  ±a  ,  (62) 

since  t+  *  0  and  x  *  ±a  [Eq.  (60)]  for  t  =  0,  and,  hence,  f+  *  f+(0)  *  0  [Eq.  (58)] 
and  v(f+)  *  v(0)  *  0. 

For  a  brief  illustration  of  the  transformations  (58)  consider  the  simple  con- 

i  ductor  motion  v(t)  ■  v  ,  x(t)  *  v  t,  t  >  0.  In  this  case,  t ,  “  (1  +  v  /c)t  and 

-  o  o  ±  o 

V 

|  t  *  t+/ (1  ^  — 2.)  f  i#e.,  f+(t+)  is  proportional  to  t  . 

i 

i 

i 

i 

! 

I 


X 


I 

ANALYTICAL  SOLUTIONS 

I 

For  the  most  general  representation  of  the  electromagnetic  fields  in  the  moving 
conductor  and  the  ambient  medium,  dimensionless  independent  and  dependent  variables 
are  introduced  by 

5  *  x/a,  t  -  Kt/a2,  C(t)  -  x(t)/a,  a  *  k  a,  v(r)  *  v(t)/v  (63) 

n  n  o 

8(c,t)  =  B(x,t)/Bo,  E(£,t)  »  E(x,t)/0cBo/a),  j(5,x)  -  j  (x,t)/ (BQ/Wa) , 

-  Y(x,t)/Bo  .  (64) 


ot  ctgot  3  "  ^  y  n  ®  1)2|3)  •  < 
n  n 


K  *  (R2  +  a2)sina  /a2[(R2  +  a2)  +  R] 
n  n  n  n  n 


M  »  uaav 


R  -  paac 


M  is  known  as  the  magnetic  Reynolds  number  of  the  conductor  with  characteristic  speed 

Eqs-  (65)  -  (67)  indicate  that  M/(L  +  R)  determines  the  order  of  the  ratio  B/Bq 
of  induced  and  external  magnetic  fields.  The  steady-state  induction  in  moving  con¬ 
ductors  -^*12)  determined  only  by  M. 

-1/2 

R  is  a  new  dimensionless  group  which  involves  the  velocity  of  light  c  *  (ep) 
so  that  R  >>  M  for  vq  «  c.  R  has  the  physical  meaning  of  a  "magnetic  Reynolds  num¬ 
ber"  of  the  ambient  nonconducting  space,  i.e.,  R  is  a  coupling  parameter  between  the 
conductor  (0  <  o  <  »)  and  its  external  medium  (a  *  0),  which  determines  the  magnitudes 
of  the  external  electromagnetic  transients  [Eqs.  (75)  -  (76)]. 

The  net  electric  current  flowing  through  the  conductor  is  per  unit  width  An  *  1 

£(t)+1 

I(X)  -  /  ja,T)d£  -  Mv(t)/(1  +  R) 


•  T,/I\  (T-Tf) 

-  2M  l  sinonKn  /  dy^  e  dt*  ,  t  >  0 


.  (72) 


2 .  External  Solutions 


For  the  9pace  surrounding  the  moving  conductor,  the  wave  functions  ^+(x»t)  of 
the  electromagnetic  transients  are  by  Eqs.  (60), (63)  -  (64)  in  dimensionless  repre¬ 
sentation 


'-r  -1 


mt  +  *  *2 (A  &)  v(*±(t  + 


Mr  *ak{T  +  dv(T’)  V’  , 

*  i  "*  j  dT 


£C0  ±  1  S  c  *  ±(1  +  Rt)  ,  t  >  0,  (73) 


where 


►+{t  +  s  (K/a2)f  (t  + 


(74) 


By  Eqs.  (18)  -  (19) ,  and  (73)  the  dimensionless  solutions  for  the  electromagnetic  field 
outside  the  moving  conductor  are: 


B±(£,t)  -  1  +  *±(t  + 
-  1 


£CO  ±  1  *  C  s  ±(1  +  Rt),  t  >  0 
1  +  Rt  <  ±£  <  “  >  t>0 


(75) 


and 

E±(5,t)  -  *R*+(t  +  £-|-^)  ,  5(0  ±  1  5*  £  S  ±(1+Rt)  ,  t>0, 

-0  ,  1  +  Rt  <  ±£  <  «  ,x>0.  (76) 
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DISCUSSION 


The  magnetic  Reynolds  number  is  M  *  paavQ  ^  1  for  conductors  (depending  on  the 
parameters  a,  a,  and  v^) ,  whereas  the  coupling  number  R  *  p^aac  >>  1  is  large  for  macro¬ 
scopic  (a)  conductors.  In  the  most  general  case,  electromagnetic  induction  in  a  mo¬ 
ving  conductor  is  determined  both  by  U  and  R,  where  M  «  R  since  vq  <<  c. 

4  -1  -4 

1.  Case  R  »  1.  Por  most  macroscopic  conductors,  it  is  a  >  10  0  /m  and  a  >  10  m 

—7  8  2 

(y  *  4ttx10  Vsec/Am,  c  *  3x10  m/sec)  so  that  R  >  10  ,  and  by  Eqs.  (68)  and  (69) 

aQ  5  mr  ,  R  »  n*  ,  n  ■  1,2,3,  ...  ,  (77) 

Kn  *  ("1)n+1(nir>’lR~1  *  R  »  nir  ♦  n  -1,2,3,...  ,  (78) 

It  is  seen  that  the  "velocity"  and  "acceleration"  fields  in  Eqs.  (65)  ~  (67)  are  of  the 
same  order  since 


M/(l  +  R)  -  M (+K^)  -  M/R  =  vQ/c  «  1  ,  R  »  1  .  (79) 

In  the  magnetic  field  solution  (65),  the  dominant  term  is  the  external  Bq  *  1  >>  M/R, 
and  in  the  electric  field  solution  (67)  the  dominant  term  is  the  motion  induced  field 
|Mv80|  >>  |/|  >  if  R  >>  1*  Although  significant  electric  fields  are  Induced  in  the  con¬ 
ductor,  the  induced  magnetic  field  is  small,  B  ~  M/(l  +  R)  <<  B  ■  1  if  R  »  1,  and, 
therefore,  the  current  density  /  ■  38/35  is  small,  too. 

2.  Case  R  »  In  an  actual  experiment,  R  -  ®  can  never  be  reached  but  only 
asymptotically  approached.  In  this  hypothetical  situation,  the  conductor  (Eqs.  (65)  - 
(67)]  and  external  (Eqs.  (75)  -  (76)]  solutions  reduce  to: 


and 
8(5, t) 


8te,t)  -  1  ,  jOi ,T)  -  0,  E(5,t)  -  4<\>(t),  t(T)  -  1  <  l  <  c(x)  +  1,  T  >  0, 

(80) 

1  +  OfMR*1],  E±  (£ ,  t  )  -  MR0  »»  0  ,  C  (t)  ±1^5*  ±(1+Rt)  ,  t  >  0 
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* 


(81) 


since 


lima  *  nrr  ,  lim  K  *  0  ,  lim  RK  *  (-l)n+1/n7T ,  n  *  1,2,3,...  .  (82 

R-XO  n  Rxo  n  R-x°  n 

It  should  be  noted  that  for  R  «,  only  B+  =  0  but  E+  4  0,  i.e.,  the  external 
electric  transients  are  (behind  their  wave  fronts)  of  the  same  order  of  magnitu- 
de  as  the  electric  field  E  in  the  conductor  .  Since  Vxg^  *  c  3E+/3t,  the  elec¬ 
tric  transients  E+  always  coexist  with  (no  matter  how  small)  magnetic  transients 
B+ .  For  this  reason,  the  limit  R  *  00  has  no  physical  meaning,  quite  apart  from 
the  fact  that  always  R  <  «  for  a,a,c  <  00 .  The  conventional  boundary  condi tions^’^*^ 
for  electromagnetic  diffusion  processes  in  conductors  Imply  R  *  ®  and  E+  =  0,  and 
are,  therefore,  physically  incorrect. 

As  an  explanation  it  Is  noted  that.  In  conductors,  magnetic  field  diffusion 
is  a  nonrelativistic  process,  as  is  electric  conduction,  j  *  cr(E  +  vxB).  The  elec¬ 
tric  transients  E+  in  vacuum  must  be  of  the  same  order  as  the  electric  field  E  in  the 
conductors,  E+  ~  E  -  MR^,  since  otherwise  the  tangential  electric  field  would  not  be 
continuous  across  the  conductor  surface.  Cta  the  other  hand,  the  external  magnetic 
transients  are  small  of  order  MR  ^  *  U(\i^oa)  ^c  ^  since  the  magnetic  field  energy 
flows  with  the  speed  of  light  in  the  ambient  space.  The  deeper  physical  reason  for 

the  external  electromagnetic  transients  is  to  be  seen  in  the  conservation  laws  for 

9) 

electromagnetic  energy  and  momentum,  which  follow  from  Maxwell's  equations 

3.  Diffusion  Approximation.  It  is  known  that  Maxwell's  equations  with  dis- 

placement  current  and  the  nonrelativistic  Ohm's  law,  j  «  c(E  +  v*B) ,  combine  to  a 

9) 

hyperbolic  diffusion  equation  for  the  magnetic  field  B  in  conductors  ,  which  reads 
in  the  considered  one-dimensional  case  with  r-independent  conductor  velocity  v(t) 


32B  _1_  SB  v_  3B 
3t2  TR3t  TR  3x 


2  32B 


3x 


90 


■wy-'URgy 


(83) 


where 


X  =  e  /o  <<  1  (84) 

K.  O 

is  the  field  relaxation  time,  which  is  extremely  small  for  conductors.  Eq.  (83) 

reduces  to  the  parabolic  diffusion  equation  in  the  limit  t  <<  1: 

K 


3B  3B 
3t  ^ 


2 

2  a  B 

T«c  17 


te  «  1 


(85) 


The  parabolic  diffusion  equation  is  an  excellent  approximation,  since  the  relaxa¬ 
tion  time  of  conductors  is  very  small,  <<  !•  By  Eqs.  (84)  and  (85),  the  field 
relaxation  time  xR  and  the  diffusion  time  xD  are  interrelated  by 


-1 


2/  2 
Trc  /a 


u  aa 
o 


(86) 


where  a  is  the  extension  of  the  conductor.  For  conductors,  the  diffusion  time  is 
relatively  large  if  a  is  not  too  small,  i.e.,  x  >>  x  . 

“  D  K 

2  2 

Comparison  of  the  neglected  term  9  B/Dt  with  the  leading  (0  <  |v|  <<  c) 
second  and  fourth  terms  of  Eq.  (83)  reveals  the  relation  of  the  parabolic  diffusion 
approximation  to  the  new  coupling  number  R  =*  \iqo ac: 


i32B|  ,|  -1  3B,  |  3  2B.  2  32Bi  TR 

y ' ,|T*  **• '  17  1  7?  '  “ 


eoyo 


n~2 


D  (uQaa ) 


(87) 


This  result  again  confirms  the  validity  of  the  parabolic  diffusion  eauation  for 

conductors,  for  which  R  ■  uQoac  >>  1.  More  important,  Eq.  (87)  demonstrates  that 

2  2  -2 

the  neglected  relativistic  term  3  B/3t  is  small  of  order  R  <«<  1,  whereas  the 
calculated  electromagnetic  fields  in  the  conductor  are  of  order  B  ~  E  ~  MR0  [Eqs. 
(65)  -  (67)],  and  the  external  electromagnetic  transients  are  of  order  B  *  MR  ^ 


and  E+  -  MRU  [Eqs.  (73)  -  (76)], 
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Thus,  consistent  electromagnetic  field  solutions  for  the  regions  in  and 
outside  of  the  accelerated  conductor  have  been  obtained  within  the  parabolic 
diffusion  approximation,  which  satisfy  the  boundary  conditions  for  the  conti¬ 
nuity  of  the  tangential  electric  and  magnetic  fields.  If  an  accuracy  of  higher 
-2 

order  R  is  to  be  achieved,  then  the  hyperbolic  diffusion  equation  (83)  has  to 
be  used.  For  non-relativistic  conductor  motions,  )v)  <<  c,  however,  an  accuracy 
of  order  R  ^  is  completely  sufficient.  The  mathematical  advantages  of  the  para¬ 
bolic  diffusion  equation  become  obvious  if  it  is  used  in  connection  with  (time- 
dependent)  moving  boundary  conditions  (magnetic  flux  compressors,  electromagnetic 
induction  generators,  etc.),  which  are  extremely  difficult  to  treat  for  the  hyper¬ 
bolic  diffusion  equation. 
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ABSTRACT 


The  heating  rate  of  electrons  in  laser-irradiated  plasmas  is  derived 
from  the  quantum-mechanic ally  extended  Vlasov  equation.  The  heating  of 
electrons  is  shown  to  be  due  to  (i)  the  electron  velocity  space  diffusion 
by  the  plasma  waves,  which  are  induced  by  the  laser  field  (classical  tur¬ 
bulent  heating),  (ii)  the  microscopic  quantum-mechanical  interaction 
between  each  electron  and  the  laser  field  in  the  electric  field  of  the 
plasma  waves  (collisionless  inverse  bremsstrahlung) ,  and  (iii)  collisio- 
nal  inverse  bremsstrahlung  in  the  static  plasma  field.  A  quasilinear  equa¬ 
tion  for  the  changes  of  the  electron  distribution  function  by  these  three 
concomitant  mechanisms  is  derived  and  the  heating  rates  of  the  electrons 
are  calculated.  It  is  shown  that  the  collisionless  inverse  bremsstrahlung 
is  dominant  for  hot  plasmas  in  most  experimental  cases.  The  so-called  anomalous 
heating  is  identified  as  being  due  to  collisionless  inverse  bremsstrahlung, 
and  is  explicitly  expressed  in  terms  of  the  wave  number  of  the  plasma  wave 
and  the  plasma  parameters,  when  the  turbulence  is  due  to  the  ion  acoustic 
instability  and  is  stabilized  by  ion  trapping. 


94 


AD-A112  099  NAVAL  WEAPONS  CENTER  CHINA  LAKE  C *  F/G  ?0/3 

EXPLOSION  DRIVEN  MAGNCTOOASOTNAMK  FLOWS  WITH  HIGH  MAGNETIC  RFT— ETC(H) 
DEC  61  HE  WILHELM 


UNCLASSIFIED 


INTRODUCTION 


If  a  laser  beam  is  injected  into  plasma,  the  energy  of  the  laser  field 
is  transferred  to  the  electrons  by  the  following  mechanisms: 

(i)  In  part,  the  laser  radiation  is  transformed  into  plasma  waves  through 
interactions  with  electrons  or  the  already  existing  plasma  waves ,  Then, 
energy  flows  from  plasma  waves  to  electrons  by  Landau  damping,  which  diffuses 
the  electrons  in  velocity  space.  This  mechanism  is  known  not  only  as  the  only 
one  heating  mechanism  caused  by  the  plasma  waves  but  also  as  the  responsible 
mechanism  for  the  anomalously  large  heating  of  collisionless  plasmas  by  la¬ 
sers  [1].  Hence,  this  mechanism  has  been  called  either  turbulent  heating 

or  anomalous  heating  by  other  authors  [1,2,3].  However,  we  call  it  "classi¬ 
cal  turbulent  heating"  since  this  mechanism  turnes  out  to  be  neither  the  sole 
heating  mechanism  of  the  plasma  waves  nor  the  dominant  cause  for  the  so-called 
anomalous  heating. 

(ii)  The  direct  absorption  of  the  laser  radiation  by  the  individual  electrons 
under  the  influence  of  the  plasma  waves  (collisionless  inverse  bremsstrahlung) 
[4]. 

(iii)  The  same  process  as  collisionless  inverse  bremsstrahlung,  but  under  the 
influence  of  the  static  plasma  field  (collisional  inverse  bremsstrahlung)  [5]. 

These  three  mechanisms  are  independent  but  can  occur  concomitantly  as 
one  readily  sees  from  the  quantum-mechanically  extended  Vlasov  equation  (Kim  - 
Wilhelm  equation) . 

The  first  systematic  theory  of  the  dissipation  of  an  electromagnetic 
wave  in  plasmas  has  been  given  by  Dawson  and  Oberman  [6].  However,  since  their 
theory  is  based  on  the  hydrodynamic  momentum  equation  (force  equation)  of  the 
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electrons  only,  the  heating  of  the  electrons  cannot  be  treated  by  their  theory 
in  principle  (only  the  over-all  ohmic  dissipation  of  the  electromagnetic  wave 
energy,  including  scattering,  can  be  determined).  Furthermore,  the  application 
of  their  theory  is  limited  in  the  weak  laser-intensity  regime  where  anomalous 
heating  by  laser-induced  plasma  waves  does  not  exist. 

Recently,  Manheimer  [7]  has  conjectured  an  equation  for  the  heating  of 
the  electrons  in  the  laser-irradiated  plasma  by  adding  a  phenomenological  term 
for  the  anomalous  heating  to  the  energy  equation  of  the  electrons  [7,8]*  The 
latter  appears  to  have  been  guessed  from  the  Vlasov  equation  without  taking  in 
account  properly  the  laser  field.  The  term  for  the  so-called  anomalous  heating, 
which  is  supposed  to  consider  heating  by  collisionless  inverse  bremsstrahlung, 
is  purely  phenomenological,  i.e.,  is  not  derived  from  first  principles. 

For  these  reasons,  herein  the  equation  for  the  heating  rate  of  the  elec¬ 
trons  of  the  plasma  in  the  field  of  high-frequency  electromagnetic  waves  is 
derived  from  the  quantum-mechanically  extended  Vlasov  equation.  The  quantum- 
mechanically  extended  Vlasov  equation  is  obtained  from  the  ordinary  classical 
Vlasov  equation  by  replacing  the  term  for  the  change  of  the  electron  distri¬ 
bution  due  to  acceleration  (or  deceleration)  of  the  electrons  by  the  laser 
radiation  through  a  term  representing  the  classical  limit  of  the  change  in 
the  electron  distribution  due  to  absorption  (or  emission)  of  laser  quanta  [9]. 
The  equation  for  the  heating  rate  of  electrons  here  gives  the  detailed  features 
as  to  how  and  what  amount  of  the  energy  of  the  laser  radiation  is  transferred  to 
the  electrons.  Especially  the  dominant  mechanism  for  the  anomalous  heating  of 
collisionless  plasmas  by  lasers  is  clearly  revealed  by  this  equation. 
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BASIC  EQUATIONS 


We  begin  by  adopting  the  get  of  Vlasov  equations  to  describe  the  change 
of  the  electron  distribution  in  the  (r,v)  space,  and  the  time-dependent  Schroe- 
dinger  equation  to  describe  the  dynamics  of  the  individual  electrons  in  the  r- 
space.  The  ions  are  regarded  as  a  set  of  randomly  distributed  fixed  scatters. 

We  assume  that  the  laser  light  is  a  circularly  polarized  electromagnetic 
wave  propagating  in  the  z-direction.  The  spatial  dependence  of  this  electro¬ 
magnetic  wave  is  neglected  since  the  wave  length  of  the  electromagnetic  wave 
is  assumed  to  be  much  longer  than  that  of  the  plasma  waves.  Under  these 
assumptions,  the  field  equations  of  the  plasma  are  (Gaussian  units): 


3 1 


+  v* 


11 

-+ 

3r 


-r  (t 


m 


e.m. 


- 

3v 


4ir[e/f  d  v  -  Ze£fi(r  -  r.)] 

j  J 


£ti  [(l/2m)(|  V  +  |I)2  -  e*]V 

where 

E  -  .Ill 

e.m.  c9t 


,  (i) 

,  (2) 

,  (3) 

,  (4) 


A(t)  ■  Aq(x  co suit  +  y  slnut)  .  (5) 

$(r,t)  is  the  potential  of  the  plasma  field  (the  plasma  waves  plus  the  static 
plasma  field).  The  subscript  e.m.  is  an  abbreviation  for  "electromagnetic", 
and  the  other  symbols  have  the  usual  meaning. 
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The  potential  of  the  plasma  field  can  be  expanded  in  a  double  Fourier 


series, 

Hr.o  -  l  .  (6) 

la 

(For  the  case  that  the  static  plasma  field  is  negligible,  <j>(r,t)  can  be  ex¬ 
panded  in  a  single  Fourier  series  since  in  this  case  1c  and  ft  are  interrelated 
by  the  dispersion  relation).  Since  <t>(r,t)  is  reali  we  have  <j>(-1c,-ft)  *  $  (1c,  ft)*. 

After  substituting  Eqs.  (5)  and  (6)  into  Eq.  (3),  we  solve  for  the  tran¬ 
sition  probability  amplitude  from  the  resulting  equation  by  means  of  first- 
order  perturbation  theory  [4,5]  where  <J>  is  the  pertubing  potential.  The  tran¬ 
sition  probability  amplitude  from  a  state  1  with  momentum  p^  to  a  state  2  with 
momentum  p^  of  the  electron  is  [4,5] 


a(l  -  2)  -  l 


ein6Jn(X/o.) 


sin[(^  -  ft  -  nw)T] 
(e/tf)  -  ft  -  no> 


where 


eE  k,  /mu) 
o  i 


(p2  -  P^)/2m 


6  is  the  azimuthal  angle  of  k,  i.e.,  k  *  (k^,6,k  )t  and  c*  is  a  real  constant, 
which  drops  out  since  exp(-icf)  does  not  affect  the  transition  probability. 

The  subscript  1  refers  to  the  direction  perpendicular  to  the  direction  of  pro¬ 
pagation  of  the  electromagnetic  wave,  and  Eq  is  the  electric  field  amplitude 
of  the  electromagnetic  wave. 

From  Eq.  (7)  we  see  that  the  transition  probability  per  unit  time  for  the 
transition  from  state  1  to  state  2  (with  the  absorption  or  the  emission  of  pho¬ 
tons)  is 
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(9) 


T (1  -  2)  =■  T (v1»v2>  -  |a(l  -  2) |  /2T 
2  00 

-  I  J^(X/u)6(e  -  m  -  ntlw) 

£  Q  n=-« 

*  n^O 

In  the  derivation  of  Eq.  (9),  we  assumed  that  there  are  no  phase  relations  be¬ 
tween  different  components  of  the  plasma  field. 

The  change  of  the  electron  distribution  due  to  the  absorption  or  emission 
of  photons  of  the  electromagnetic  wave  [4,5]  is 

[3f(v2)/3t]e  m  -  £[7(1  -  2)f(vx)  -  7(2  ^  l)f(v2>]  =  [f  (^)  -  H^)} 


where 


V-  -  p-/m  *  2Tfftm  (n  /L  ,n  /L  ,n  /L  )  ,  n  ,n  ,n  *  integers, 

11  x x*  y  y*  z  z  x  y9  z  ° 


is  quantized  (Lx,L^,Lz  are  the  x-,y-,z-dimensions  of  the  plasma,  respectively). 

In  Eq.  (1),  (e/ra)^  *3f/3v  is  the  rate  of  change  in  the  number  of  electrons 

e .  m. 

per  unit  volume  in  the  (r,v) space  due  to  acceleration  (or  deceleration)  of  elec¬ 
trons  by  the  laser  field.  From  the  quantum-mechanical  viewpoint,  the  acceler¬ 
ation  of  an  electron  by  the  spatially  uniform  laser  field  is  possible  only  by 
absorption  (or  emission)  of  the  laser  photons  in  the  presence  of  another  field 

which  is  the  plasma  field  in  our  situation.  Hence,  (e/m)l  *3f/3v  is  equal 

e  .m 

to  the  classical  limit  (h  -►  0)  of  the  rate  of  change  in  the  mumber  of  electrons 
per  unit  volume  in  the  (r,v)  space  due  to  the  absorption  (or  emission)  of  the 
laser  photons  by  the  electrons  under  the  influence  of  the  plasma  field  ("inverse 
bremsstrahlung") .  Accordingly,  by  using  Eqs.  (1),  (3)  -  (6),  and  the  same  meth¬ 
ods  as  in  Ref.  5  and  6,  we  obtain  the  quantum-mechanical  Vlasov  equation,  which 
is  valid  up  to  second  order  in  : 
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(11) 


3f  ,  t,  31  .  e  _A  3f  ,3f. 

sT  + v'  (i7>i.b. 

9r  3v 


where 


b  =  lim  I  T^'  -’•  v)[f(v')  -  f(v)] 

0  ^(*1  (\j  fi  I IKU  \  n 

I  !!**.<»  1 2  l  k 


^  *•  W  a  v,  a  v.  ■  t 

m  ■£  Q  n*-<*  tc  1c  K 

n^O 


k  2 
3v* 


(12) 


lung”.  In  the  derivation  of  Eq.  (12),  v  +  *  -v+,  and  <j>(-k,-a)  *  <j>(1c,fl)*  have 

been  used  to  cancel  out  the  terms,  which  change  their  signs  as  (1c,ft,n)  changes 

to  (-1c,-ft,-n),  in  the  summation  over  1c, ft,  and  n. 

Equation  (11)  is  the  desired  Vlasov  equation  of  the  electrons  in  the 

electromagnetic  field  and  the  plasma  field,  in  which  the  term  (9f/3t).  ,  con- 

i.b. 

siders  all  interactions  of  the  electrons  with  electromagnetic  field.  It  should 
be  noted  that  the  presented  derivation  of  Eq.  (11)  is  based  on  first  principles 
(Vlasov  and  Schroedinger  equations  in  electromagnetic  and  plasma  fields). 

Since  Eq.  (11)  is  valid  up  to  the  second  order  in  <J> ,  we  can  apply  to  this 
equation  an  analogy  to  the  derivation  of  the  quasilinear  equation  for  the  electron 
distribution  in  the  absence  of  the  electromagnetic  wave  [10,11].  The  result, 
i.e.,  the  quasilinear  equation  for  changes  in  the  electron  distribution  in  the 
presence  of  the  electromagnetic  wave,  is  by  Eq.  (11) 


8f 

at 


(— ) 

v3t'c.t. 


(— ) 

'‘at'i.b. 


(13) 


where 


(— ) 

v3t'c.t. 
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3v  K-V  -  G  -  ly 


(14) 
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is  the  change  in  the  electron  distribution  associated  with  the  diffusion  in 
the  electron-velocity  space  due  to  the  acceleration  of  the  electrons  by  the 
electric  field  of  the  plasma  waves  (or  the  change  in  the  electron  distribu¬ 
tion  by  classical  turbulent  heating).  The  subscript  c.t.  is  an  abbreviation 
for  "classical  turbulent  heating",  and  y  >  0  is  infinitesimal.  In  Eqs.  (13) 
and  (14) ,  the  spatial  average  of  the  actual  distribution  function  is  desig¬ 
nated  as  f,  and  the  same  simplification  in  notation  is  adopted  in  the  fol¬ 
lowing,  where  there  is  no  chance  of  confusion. 
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HEATING  RATE  OF  ELECTRONS  IN  COLLISIONLESS  PLASMA 

The  heating  rate  of  the  electrons  per  unit  volume  under  the  influence 
of  the  external  electromagnetic  and  internal  plasma  fields  is  by  microscopic 
definition 


3<e> 

8t 


1 "^\2r/3f\  ,  /3f\  i  /3<e>x 

d  v-jNb(v  -  u)  K^)1-b-  +  <3t-)c.t.l  ' 


3t  c.t. 


(15) 


where  N  is  the  unperturbed  electron  density.  In  Eq.  (15)*  the  first  term  is 
the  heating  rate  of  the  electrons  per  unit  volume  by  inverse  bremsstrahlung* 
and  the  second  term  is  that  by  classical  turbulent  heating*  i.e.,  the  rate  of 
energy  transfer  to  the  electrons  from  the  plasma  waves  by  velocity-space  dif¬ 
fusion. 

We  model  the  spatial  average  of  the  electron  distribution  as 


f(v>  -  (2,»y3,2OT(-^%l  jjl  ♦  UHkSlf  (l^^l2  .  3)}  (16) 

2v  *  3Nrav  v  Q4  7 

e 


where  the  hydrodynamic  velocity  u  and  heat  flux  y  are  induced  by  the  laser  field. 
In  the  appendix  we  calculate  the  second  term  of  Eq.  (15)  as 


(^)t ,  -  <f>1/2  *4 1  jMdiili  Mpi-w-tVia  -  5-S) 
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(17) 


The  calculation  of  the  first  term  by  means  of  Eqs.  (12),  (15),  and  (16) 
is  elementary,  and  gives 
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Equations  (17)  and  (18)  can  be  used  irrespective  of  the  origin  of  the 
plasma  field  as  long  as  the  strength  of  this  field  is  so  small  that  the  resul¬ 
ting  perturbation  of  the  electron  distribution  can  be  adequately  treated  by 
means  of  first-order  perturbation  theory.  Bearing  this  in  mind,  we  reduce 
the  expression  of  the  heating  rate  to  a  simpler  form  for  the  cases  of  small 
and  large  X/co. 

(i)  X/o)  <<  1:  In  this  case,  the  laser- induced  current  and  heat  flux 
can  be  neglected.  Then,  by  means  of  the  same  approach  used  to  derive  Eqs.  (1) 

-  (15)  of  Ref.  6,  we  ontain  from  Eqs.  (1)  and  (2) 


.  ,■*  x  (2ir)^  v  i£*rr^/*x  -i5  iwt  .  i6  -iwt, 
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(19) 

(20) 

(21) 

(22) 
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Here,  D(k,w)  is  the  longitudinal  dielectric  constant,  V  is  the  volume  of  the 


plasma  system,  83  Wp/Ve  is  the  Debye  wave  number,  the  plasma  frequency, 
and  W(s)  is  the  plasma  function  given  in  the  appendix. 

Designating  by  $ ^ (r)  the  Coulomb  field  potential  at  r  due  to  an  ion  at 
Ty  and  by  XD  the  Debye  length*  the  first  term  of  Eq.  (19)  becomes 


$(r)  -  Z*j(r)  -  I  ■  Z*  exp(-|r-r  |/X  )  ,  (23) 

Hence,  the  first  term  of  Eq.  (19)  is  the  static  Coulomb  field  potential,  re¬ 
duced  by  the  usual  Debye  shielding,  acting  on  an  electron  due  to  all  ions  in 
the  plasma.  This  static  plasma  field  gives  rise  to  collisional  inverse  brems- 
st rah lung  but  does  not  contribute  to  classical  turbulent  heating.  The  second 
and  third  terms  are  the  potential  of  the  laser-induced  plasma  waves,  which  can 
contribute  both  to  collisionless  inverse  brensstrahlung  and  classical  turbulent 
heating. 

From  Eqs.  (18)  -  (21)  and 


Jn(1) 
n  a) 


(J_)2lnl  for  x/m  «  x 

(n!)2  2w 


it  is  seen  that  the  collisionless  inverse  bremsstrahlung  is  of  the  fourth 
order  in  X/o>,  and  can  be  neglected  compared  with  the  collisional  inverse  brems¬ 
strahlung,  which  is  of  the  second  order.  Hence,  up  to  the  lowest  order  in  x/w, 
we  obtain  from  Eqs.  (18)  -  (21),  and  (24), 
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.  (25) 


104 


The  expression  for  <)>j(r)  given  in  Eq.  (23)  shows  that  <j>  (r)  decreases 
very  fast  to  zero  as  |^-r_.|/XD^  This  decrease  is  mathematically  brought 
about  by  representing  discrete  phenomena  through  smoothly  continuous  functions; 
<j>j(r)  is  in  reality  zero  for  [  r-r^  |  >>  This  corresponds  to  S (k)  *  0  for 

k  <<  k^.  Considering  this  behavior  of  S(k)  and  Eq.  (20),  we  can  explicitly 


evaluate  Eq.  (25)  as 
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Where 


(k  +  k£)‘ 


U(z)  »  /  dx[x^exp(-l/2x^)/Cx^+z^)^]  ,  z  =  w  /w 

0  P 


with  ojpthe  plasma  frequency.  The  function  U(z)  is  represented  in  Fig.  1. 

In  the  derivation  of  Eq.  (26),  we  assumed  the  ion  distribution  to  be 
random  so  that 


l  exp(-ilc*r  )  |2  *  N  (28) 

J  J 

since  the  large-amplitude  ion  wave  cannot  be  induced  by  weak  laser  fields  for 
which  X/o)  <<  1  [12]. 

From  Eqs.  (17)  and  (19),  we  obtain  for  the  heating  rate  per  unit  volume 
by  classical  turbulent  heating,  up  to  the  lowest  order  in  X/u>, 


,3<e>>  -  „>  >7/2  e2N  rd\|c(t)|2  2  ,  2,„_.2,_2> 

(ir kt.  (2l°  ~  Jk  V  "  exp(_W  /2vek  > 


Eq.  (29)  can  be  expressed  more  explicitly  for  the  cases  of  hot  and  cold  plas¬ 
mas.  In  this  connection,  it  is  noted  that  electromagnetic  waves  propagate 
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Fig.  Is  U(z)  versus  z  by  Eq.  (27) 


I 

I 


in  plasmas  only  when  their  frequency  is  greater  than  the  plasma  frequency. 
Accordingly, 


v  <  U)V  /w  =  m/k_ 
e  e  p  D 


(30) 


For  hot  plasmas  for  which  >>  oj/k,  Eq.  (22)  becomes,  to  the  lowest  or¬ 
der  in  w/kv  ,  such  that  [1/D(k,0)  -  1/D(k,u))]  is  not  zero. 


D(k,u>)  =  1  +  kj/k2  -  i(*/2)1/2(k£/k2)~ 

e 


(31) 


Also,  from  >>  cj/k  and  Eq.  (30),  we  have  k  >>  k^.  This  means  that  the 
plasma  behaves,  as  far  as  the  plasma  waves  are  concerned,  like  a  system  of  free 
individual  particles  [13]. 

Substituting  Eqs.  (21)  and  (31)  into  Eq.  (29)  and  using  k  >>  k^  yields 
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By  comparison  of  Eq.  (32)  with  Eq.  (26),  it  is  recongized  that  the  hea¬ 
ting  rate  of  the  electrons  by  inverse  bremsstrahlung  is  slightly  greater  than 
that  by  classical  turbulent  heating  for  o>  » 

For  cold  plasmas  for  which  v^  <  w/k^  <<  oi/k,  Eq.  (22)  becomes,  to  the 
lowest  order  in  v^k/w  such  that  [1/D(k,0)  -  1/D(k,u>)]  in  Eq.  (21)  is  not  zero. 


D(k,o>)  »  1  -  u>^/ u>^ 
P 


Accordingly,  Eq.  (29)  reduces  to 
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(33) 


(34) 
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where 


V(z)  *  /  dx[exp(-l/2x2)/x(x2+z2)2]  .  (35) 

0 

The  function  V(z)  is  represented  in  Fig.  2. 

Equation  (34)  appears  to  indicate  that  classical  turbulent  heating  by 
the  plasma  wave  becomes  resonant  as  ui  -►  u)p.  However,  Eq.  (34)  should  be  ques¬ 
tioned  since  Eq.  (38)  is  derived  by  using  Eqs.  (19)  and  (21),  which  are  derived 
from  the  collisionless  Boltzmann  equation  (or  Vlasov  equation) ,  even  though  the 
considered  plasmas  are  cold,  and  consequently,  collisional. 

Comparison  of  the  expression  for  the  heating  rate  by  inverse  bremsstrah- 
lung,  i.e.,  Eq.  (26)  with  that  found  by  Seely  and  Harris  [Eq.  (15)  of  Ref.  5] 
shows  that  our  result  is  the  same  as  theirs  except  for  a  factor  U(z)/lnA,  where 
InA  is  the  Coulomb  logarithm. 


InA  -  ln(127rNX^) 


(36) 


Since  the  Coulomb  logarithm  for  laboratory  plasmas  is  about  10  to  20  and  U(z) 

is  less  than  0.04,  the  value  of  the  heating  rate  given  by  Eq.  (26)  is  less  than 

about  four-thousandth  to  two- thousandth  of  that  given  by  Eq.  (15)  of  Ref.  5. 

This  difference  is  due  to  the  fact  that  Seely  and  Harris  do  not  consider  Debye 

shielding,  whereas  our  result  contains  Debye  shielding  brought  about  by  the 

collective  dynamics  of  the  electrons  [Eq.  (19)]* 

2 

(ii)  A/w  ■  eE^/mw  >>  1:  For  this  case,  we  can  use  no  longer  Eqs.  (19) 
-  (21)  for  the  plasma  field  potential  $  even  if  there  is  no  other  perturbation 
except  the  electromagnetic  field.  However,  we  can  readily  surmise  that  the 
collisionless  inverse  bremsstrahlung  is  stronger  than  the  collisional  one  for 
y  =  X/u>  >>  1.  This  is  because  the  second  and  third  terms  of  Eq.  (19),  which 
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give  rise  to  collisionless  inverse  bremsstrahlung,  become  important  compared 
to  the  first  term,  which  produces  collisional  inverse  bremsstrahlung,  as  y 
increases  (but  still  y  <  1) .  For  this  reason,  the  result  obtained  by  Seely 
and  Harris  for  the  effective  collision  frequency  in  the  case  y  >>  1  (Eq.  (22) 
of  Ref.  5]  is  questionable  since  they  take  into  account  only  collisional  in¬ 
verse  bremsstrahlung. 

Except  for  increasing  y  (usually  by  increasing  the  intensity  of  the  elec¬ 
tromagnetic  field),  the  collisionless  inverse  bremsstrahlung  becomes  more  impor¬ 
tant  compared  with  the  collisional  one  with  increasing  electron  temperature. 

This  is  because  the  mean  free  path  of  the  electron  becomes  larger  with  increasing 
electron  energy.  Accordingly,  the  collective  effects  which  are  enhanced  by  a 
large  number  of  instabilities  (not  caused  by  the  local  electromagnetic  field) 
become  more  dominant  than  the  effect  caused  by  the  static  ions.  Hence,  even  if 
y  <  1,  the  collisionless  inverse  bremsstrahlung  may  dominate  the  collisional  one 
for  very  hot  plasmas  if  there  are  causes  of  turbulence  other  than  the  local  elec¬ 
tromagnetic  field.  The  latter  case  will  be  treated  in  the  next  section. 

For  y  *  X/w  >>  1,  we  can  write  approximately  [51 


„LJ»6)S  (vt '  H584  •  i'5  (vt  -  H** +  s  <v£  -  1  •  <37> 

n^O 

Substituting  Eq.  (37)  into  the  first  relation  of  Eq.  (18),  we  find  for  hot 
plasmas,  >>  ft/k,  X/k, 
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Comparison  of  Eq.  (38)  with  Eq.  (17)  indicates  that  the  heating  of  the 
electrons  by  collisionless  inverse  bremsstrahlung  is  greater  than  that  by 
classical  turbulent  heating  where  the  following  conditions  hold:  u  >>  1  and 
ve  **  */k  (hot  plasmas). 
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ANOMALOUS  HEATING  RATE 


Heating  by  collisionless  inverse  bremsstrahlung,  which  is  the  main  heating 
mechanism  in  laser  fusion,  is  usually  called  "anomalous"  heating.  This  phenomenon 
requires  a  correct  quantum-mechanical  treatment  of  collisionless  inverse  brems- 
st rah lung,  as  distinguished  from  the  deficient  classical  approaches  in  Refs.  [6,7, 
8,12,13].  From  the  experimental  point  of  view,  the  heating  rate  is  "anomalous*'  in 
collisionless  plasmas  since  it  is  unexpectively  large  compared  with  the  classical 
theories,  which  consider  collision-dominated  phenomena.  In  the  following,  the 
significance  of  collisionless  inverse  bremsstrahlung  is  demonstrated  by  comparison 
with  the  classical  theory. 

Equation  (38)  is  valid  irrespective  of  the  causes  of  the  electrostatic  field. 

In  laser  fusion,  it  is  believed  that  the  electrostatic  field  is  developed  by  the 

ion  acoustic  instability  driven  by  the  return  current  [7].  For  this  case,  it  is 

known  from  computer  simulations  that  the  k-spectrum  of  the  electrostatic  field  is 

flared  out  in  a  cone  of  45®  to  the  direction  of  the  return  current  [7,15-17].  Hence, 

we  can  assume  that  k^  is  of  the  order  of  k.  Then,  the  condition  y  >  1  corresponds 

12  2 

to  a  laser  intensity  of  I  >  10  watt /cm  for  a  neodymium  laser*  With  presently 
available  laser  powers,  this  condition  is  generally  satisfied  in  most  experiments 
so  that  Eq.  (38)  is  applicable  to  them. 

Both  computer  simulations  and  analytical  theories  Indicate  that  the  ion 
acoustic  instability  is  "stabilized"  by  ion  trapping  [18-22].  Hence,  its  Intensity 
level  can  be  estimated  on  the  basis  of  simple  trapping  arguments  to  be  [7] 
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where  the  weve  spectrua  of  the  turbulent  electric  field  Is  modeled  for  a  single 
wave  having  wave  number  k,  frequency  ft,  and  a  propagation  angle  of  133°  with  the 
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direction  of  the  laser  beam.  By  incorporating  the  estimate  (39)  for  the 
saturated  turbulence  level  into  Eq.  (38),  we  obtain  the  anomalous  heating 
rate  for  hot  plasmas  with  ve  >>  X/k: 
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By  Eqs.  (26),  and  (40),  the  ratio  R  of  the  anomalous  heating  rate  to  the 

heating  rate  by  the  collisional  inverse  bremsstrahlung  in  very  hot  (collision¬ 
less)  plasmas  is,  assuming  Tg  »  and  k  »  k^/2: 

R  *  NAjJ/U(z)  *  6Nd  ,  Nd  -  4irApN/3  ,  (Al) 

where  is  the  number  of  electrons  in  a  Debye  sphere  [23].  By  definition,  in 
"collisionless  plasmas"  the  collective  behavior  far  dominates  the  collisional 
effects  [23],  and  >>>  1.  Thus,  by  Eq.  (41)  we  demonstrate  the  significance 
of  our  collisionless  quantum-mechanical  theory  compared  with  the  classical  col¬ 
lisional  theories  and  explain  the  observed  anomalous  heating  rate  in  experiments. 


R  -  6N^  >>>  1. 


CONCLUSION 


Under  the  condition  that  the  field  strength  of  the  electromagnetic  wave  is 
below  the  threshold  above  which  strong  turbulence  is  developed  (so  that  the  first- 
order  perturbation  theory  would  no  longer  be  appropriate)  we  have  obtained  the  fol¬ 
lowing  results: 

(1)  From  the  set  of  Vlasov  and  Schroedinger  equations,  the  heating  rate  of  the 

electrons  in  the  electromagnetic  field  is  systematically  derived  in  terms  of  the  elec¬ 
trostatic  field  potential,  which  is  produced  by  the  static  ions  and  the  collective 
dynamics  of  electrons  in  the  self-consistent  manner.  The  Vlasov  equation  is  genera¬ 
lized  by  formulating  in  it  —1  *3f/3v  as  a  quantum-mechanical  interaction  integral. 

m  e .m. 

(2)  Especially  for  y  <  1,  the  electrostatic  field  potential  in  the  heating 
rates  is  expressed  in  terms  of  the  field  strength  Eq,  frequency  w  ,  and  wave  vector  £ 
of  the  electromagnetic  wave  for  both  hot  and  cold  plasmas. 

(3)  For  very  hot  plasmas  with  vfi»  8/k,X/k,  it  is  shown  that  the  heating  of  the 
electrons  by  inverse  bremsstrahlung  is  considerably  stronger  than  that  by  Landau 
damping  of  electrostatic  waves.  Furthermore,  it  is  shown  that  for  U  <  1  the  above 
statement  is  true  regardless  of  the  magnitude  of  the  electron  temperature. 

(4)  The  anomalous  heating  observed  in  the  experiments  is  identified  as  the 
heating  by  inverse  collisionless  bremsstrahlung,  and  explicitly  formulated  when  the 
turbulence  is  developed  by  the  ion  acoustic  instability  and  "stabilized"  by  ion 
trapping. 
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APPENDIX 


Combination  of  Eqs.  (14)  -  (16)  yields 
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Integration  over  x  in  Eq.  (Al)  is  accomplished  by  using  the  following 


formulae : 


-1/2  r  dx  x  exp (-x  /2) 


x  -  z  -  iy 


,  (A2) 


(2tt)-1/2  r~  .  zN(z) 

v  ;  x  -  z  -  iy 


(2*r1/2  /~dxxx-T-1r 


1  +  z  W(z) 


,  (A4) 


(2n)~1/2  f*°—  jLexp(-x2; 

i_  x  -  z  -  iy 


3  +  z2[l  +  z2W(z)] 


,  (A5) 


where  W(z)  is  the  plasma  function  [24]  given  by 


.1/0  2  ,  4  ,  nn+l  n+2 

i(-2>1^2z  exp(--y)  +  [1  -  z  +  -  ***  +  ^  (2n+l)!l"  +  J*for  lzl  <  1 


W(z)  - 


i(j)1/2z  exp(-^-)  -  t-7  +  "7+  •**  +  -  +  1.  for  |z|  >  1. 

z  z  z 


Substitution  of  Eqs.  (A2)  -  (A6)  into  Eq.  (Al)  and  obliteration  of  all  odd 


terms  for  both  k  and  ft  in  the  summation  over  k  and  ft  gives  Eq.  (17). 
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